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Abstract 

In this brief review we introduce the methods of quantum field theory out of equi- 
hbrium and study the non-equihbrium aspects of phase transitions. Specifically we 
critically study the picture of the "slow-roll" phase transition in the new inflationary 
models, we show that the instabilities that are the hallmark of the phase transition, that 
is the formation of correlated domains, dramatically change this picture. We analyze 
in detail the dynamics of phase separation in strongly supercooled phase transitions in 
Minkowski space. We argue that this is typically the situation in weakly coupled scalar 
theories. The effective evolution equations for the expectation value and the fluctua- 
tions of an inflaton field in a FRW cosmology are derived both in the loop expansion 
and in a self-consistent non-perturbative scheme. Finally we use these non-equilibrium 
techniques and concepts to study the influence of quantum and thermal fluctuations 
on the dynamics of a proposed mechanism for the formation of disoriented chiral con- 
densates during a rapid phase transition out of equilibrium. This last topic may prove 
to be experimentally relevant at present accelerator energies. 

To appear in the Proceedings of the '2nd. Journee Cosmologie', Observatoire de 
Paris, 2-4, June 1994- H J de Vega and N. Sanchez, Editors, World Scientific. 



1 



1 Introduction and Motivation 



Phase transitions play a fundamental role in our understanding of the interplay between 
cosmology and particle physics in extreme environments. It is widely accepted that the 
Universe that we see today is a remnant of several different phase transitions that took 
place at different energy (temperature) scales and with remarkable consequences at low 
temperatures and energies. Restricting ourselves to the domain of particle physics, there is 
the conjecture that a Grand Unified phase transition has taken place at about 10^^ Gev, in 
which a GUT symmetry breaks down to SU{3)^SU{2)<^U{1). Within the desert hypothesis, 
the next phase transition took place at about 240 Gev at which the remaining symmetry is 
broken down to SU{3) (8>f/(l)|em- Finally at a scale of about 100-200 Mev the confining and 
chiral phase transition ocurred. Heavy Ion colliders at RHIC and perhaps LHC will probe 
this energy region within the next few years, thus the understanding of the non-equilibrium 
properties of phase transitions becomes a more tangible problem. 

To put things in perspective, the temperature at the center of the sun is about 10"^ K ^ 
IkeV. 

There is also the possibility that the observed baryon asymmetry may be explained as a 
result of a phase transition at electroweak energy scales and non-equilibrium conditions [|15|. 

Phase transitions are an essential ingredient in inflationary models of the early universe [0- 
0. The importance of the description of phase transitions in extreme environments was 
recognized long time ago and efforts were devoted to their description in relativistic quantum 
field theory at finite temperature For a very thorough account of phase transitions in 
the early universe see the reviews by Br andenberger [p!0|] , Kolb and Turner |]1T| and Linde[|l^]. 

The methods used to study the equilibrium properties of phase transitions, like the order 
of the transition and critical temperature are by now well understood and widely used, in 
particular field theory at finite temperature and effective potentials ||13||. 

These methods, however, are restricted to a static description of the consequences of 
the phase transition, but can hardly be used to understand the dynamics of the processes 
involved during the phase transition. In particular, the effective potential that is widely used 
to determine the nature of a phase transition and static quantities like critical temperatures 
etc, is not the appropriate tool for the description of the dynamics. The effective potential 
corresponds to the equilibrium free energy density as a function of the order parameter. This 
is a static quantity, calculated in equilibrium, and to one loop order it is complex within the 
region of homogeneous field configurations in which V"{(f)) < 0, where V{(l)) is the classical 
potential. This was already recognized in the early treatments by Dolan and JackiwP]. 
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In statistical mechanics, this region is referred to as the "spinodal" and corresponds to a 
sequence of states which are thermodynamically unstable. 

At zero temperature, the imaginary part of the effective potential has been identified 



with the decay rate of this particular unstable state ||22||. 

The use of the static effective potential to describe the dynamics of phase transitions has 
been criticized by many authors 0, |l6l. These authors argued that phase transitions in 
typical theories will occur via the formation and growth of correlated domains inside which 
the field will relax to the value of the minimum of the equilibrium free energy fairly quickly. 
It is now believed that this picture may be correct for strongly coupled theories but is not 
accurate for weak couplings. 

The mechanism that is responsible for a typical second order phase transition from an 
initially symmetric high temperature phase is fairly well known. When the temperature 
becomes lower than the critical temperature, long-wavelength fluctuations become unstable 
and begin to grow and the field becomes correlated inside "domains" . The order parameter 
(the expectation value of the volume average of the field), remains zero all throughout the 
transition. These correlated regions grow as a function of time and eventually at long 
times (determined by non-linear and dissipative processes) the system reaches a state of 
equilibrium. Inside each correlation volume the field is at a minima of the effective potential. 

This is the process of phase separation, whether this process can be studied in equihbrium 
or not is a matter of time scales. 

Equilibrium or Non-Equilibrium: C'est la question! 

Consider the situation of an homogeneous and isotropic cosmological model in which the 
temperature of the universe decreases because of the expansion typically as 

m = nto) ^ (1) 

a[t) 

with a{t) the cosmological expansion factor. The cooling rate is thus determined by 

-f/T = H (2) 

where H stands for the Hubble constant, H{t) = a{t)/a{t). This determines the time 
scale for variation of temperature. If relaxational processes of the system are much faster 
than this scale then the system is in "quasi-thermo dynamic equilibrium" locally. Then a 
determination of whether an equilibrium description is valid requires the comparison of these 
two time scales. 

The typical collisional relaxation rate for a process at energy E in the heat bath is 
given by T{E) = t^^{E) ^ n{E)a{E)v{E) with n{E) the number density of particles with 
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this energy E, cr{E) the scattering cross section at this energy and v{E) the velocity of 
the incident beam of particles. The lowest order (Born approximation) scattering cross 
section in a typical scalar theory with quartic interaction is cr{E) ^ X^/E"^. At very high 
temperatures, T ^ m$, with m$ the mass of the field, n{T) ^ T^, the internal energy 
density is U/V ~ T^, and the average energy per particle is {E) ^ T, and v{T) ^ 1. 
Thus the typical collisional relaxation rate is r(T) ^ A^T. In an expanding universe, the 
conditions for local equilibrium will prevail provided that r(T) ^ H{t)^ with a{t) the FRW 
scale factor. If this is the case, the collisions occur very quickly compared to the expansion 
rate, and particles will equilibrate. This argument applies to the collisional relaxation of 
high frequency (short- wavelength) modes for which k ^ m$. This may be understood as 
follows. Each "external leg" in the scattering process considered, carries typical momentum 
and energy k ^ E ^ T > Tc. But in these typical theories ~ m.^/^/X. Thus for weakly 
coupled theories Tc ^ m$. 

To obtain an order of magnitude estimate, we concentrate near the phase transition at 
T ^ Tc ^ lO^^Gev, H ^ T^/Mpi ^ 10"^ T, this implies that for A > 0.01 the conditions 
for local equilibrium may prevail. However, in weakly coupled inflaton models of inflation, 
phenomenologically the coupling is bound by the spectrum of density fluctuations to be 
A ~ 10"^^ — 10~^^|T2|, |T^. Thus, in these weakly coupled theories the conditions for local 
equilibrium of high energy modes may not be achieved. One may, however, assume that 
although the scalar field is weakly self-coupled, it has strong coupling to the heat bath 
(presumably other fields in the theory) and thus remains in local thermodynamic equilibrium. 

This argument, however, applies to the collisional relaxation of short wavelength modes. 
We observe, however, that these type of arguments are not valid for the dynamics of the 
long-wavelength modes at temperatures below T^, for the following reason. 

At very high temperatures, and in local equilibrium, the system is in the disordered 
phase with ($) = and short ranged correlations, as measured by the equal time correlation 
function (properly subtracted) 



As the temperature drops near the critical temperature, and below, the phase tran- 
sition occurs. The onset of the phase transition is characterized by the instabilities of 
long-wavelength fluctuations, and the ensuing growth of correlations. The field begins to 
correlate over larger distances, and "domains" will form and grow. If the initial value of 
the order parameter is zero, it will remain zero throughout the transition. This process 



mf,t)<i>{6,t)) 



T'exp[-\r\/aT)] 



1 T 

—— ^ e(o) — 



(3) 
(4) 
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is very similar to the process of spinodal decomposition, or phase separation in condensed 



matter systems |1T8|, 0. This growth of correlations, cannot be described as a process in local 
thermodynamic equilibrium. 

These instabilities are manifest in the equilibrium free energy in the form of imaginary 
parts, and the equilibrium free energy is not a relevant quantity to study the dynamics. 

These long wavelength modes whose instabilities trigger the phase transitions have very 
slow dynamics. This is the phenomenon of critical slowing down that is observed experimen- 



tally in binary mixtures and numerically in typical simulations of phase transitions [P0|, |T 
The long wavelength fluctuations correspond to coherent collective behavior in which de- 
grees of freedom become correlated over large distances. These collective long-wavelength 
modes have extremely slow relaxation near the phase transition, and they do not have many 
available low energy decay channels. Certainly through the phase transition, high frequency, 
short wavelength modes may still remain in local equilibrium by the arguments presented 
above (if the coupling is sufficiently strong), they have many channels for decay and thus 
will maintain local equilibrium through the phase transition. 

To make this argument more quantitative, consider the situation in which the final tem- 
perature is below the critical value and early times after the transition. For small amplitude 
fiuctuations of the field, long wavelength modes "see" an inverted harmonic oscillator and 
the amplitude fiuctuations begin to grow as (see below) 

($fc(t)$_g(t)) ^ exp[2Wik)t] (5) 



W{k) = \JiJ?{T)-k^ (6) 

21 



/i'(T) = /i^(O) 



(7) 



for P < (T). 

In particular this situation, modelled with the inverted harmonic oscillators is precisely 



the situation thoroughly and clearly studied by Weinberg and Wu||22[| and Guth and Pi||23| . 

The time scales that must be compared for the dynamics of these instabilities, are now 
the growth rate T{k) ^ \J ^^^{T) — k"^ and the expansion rate H ^ T'^/Mpi ^ 10~^ T , if 
the expansion rate is comparable to the growth rate, then the long wavelength modes that 
are unstable and begin to correlate may be in local thermodynamic equilibrium through the 
cooling down process. Using Tc ~ /i(0)/-\/A, we must compare [1 - (£)2]5 to 10-Vv^. 
Clearly for weakly coupled theories, or very near the critical temperature, the growth rate 
of the unstable modes will be much slower than the rate of cooling down, and the phase 
transition will be supercooled, similarly to a "quenching process" from a high temperature. 



5 



disordered phase to a supercooled low temperature situation. For example, for A ~ 10^^^ the 
growth rate of long wavelength fluctuations, is much smaller than the expansion rate even 
for a final temperature Ty ^ 0, and the long wavelength modes will be strongly supercooled. 

The purpose of this brief review is to introduce the methods of non-equilibrium quantum 
statistical mechanics to study the dynamics of phase transitions, and apply them to relevant 
questions in cosmological phase transitions. 

Among the topics that we study here are: the "rolling" equations and the corrections 
introduced by quantum fluctuations, the process of phase separation via formation and 
growth of correlated domains, the effective equations of motion in Friedmann- Robertson- 
Walker cosmologies, including quantum and thermal corrections. We also study the proposal 
of the formation of chirally disoriented domains, which may have experimental relevance in 
P — P and heavy-ion collisions. 



2 Statistical Mechanics out of equilibrium: 

The generalization of statistical mechanics techniques to the description of non-equilibrium 
processes in quantum field theory, have been available for a long time[^]-[^| but somehow 
had not yet been accepted as an integral part of the available tools to study field theory in 
extreme environments. We thus begin by presenting a somewhat pedagogical introduction 
to the subject for the non-practioner. 

The non-equilibrium description of a system is determined by the time evolution of the 
density matrix that describes the system. This time evolution (in the Schrodinger picture) 
is determined by the quantum Liouville equation: 

in^-p^ = [H{t),p{t)] (8) 

A non-equilibrium situation arises whenever the Hamiltonian does not commute with the 
density matrix. Here we allow for an explicit time dependent Hamiltonian (for example if 
the system is in an external time dependent background). 

The formal solution of the Liouville equation is: 

p{t) = U{t,t^)p{t^)U-\t,t^) (9) 

with p{to) the density matrix at some initial time to that determines the initial condition for 
the evolution, and with U{t,to) the time evolution operator. 

In most cases of interest the initial density matrix is either thermal or a pure state 
corresponding to the ground state of some initial Hamiltonian. In both cases the initial 
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density matrix is 

p{to) = e-^^' (10) 

the ground state of Hi can be projected out by taking the hmit f3 — > oo. Non-equihbrium 
will thus arise whenever the evolution Hamiltonian does not commute with the initial Hamil- 
tonian Hi. It is convenient to introduce then a time dependent Hamiltonian H{t) such that 
H{t) = Hi for — oo < t < to and H(t) = Hevoiit) for t > to. with Hevoiit) the evolution 
Hamiltonian that determines the dynamics. This corresponds to the assumption that the 
system has been in equilibrium up to to and will evolve out of equilibrium thereafter. Real 
time equilibrium correlation functions can be obtained if H is constant in time for all times. 
The expectation value of any operator is thus 

(0)(t) = Tr [e-'^'^^U-\t)OU{t)] /Tr [e"'^'-^'] (11) 

This expression may be written in a more illuminating form by choosing an arbitrary time 
T < for which U{T) = exp[—iTHi] then we may write exp[—f3iHi] = exp[~iHi{T — i(3i — 
T)] = U(T — iPi,T). Inserting in the trace U^^{T)U{T) = 1, commuting U^^(T) with pi 
and using the composition property of the evolution operator, we may write (0) as 

{O) (t) = Tr [U{T - t(5i, t)OU{t, T)] /Tr [U{T - tf3„ T)] (12) 

The numerator of the above expression has a simple meaning: start at time T < 0, evolve 
to time t, insert the operator O and evolve backwards in time from t to T < 0, and along 
the negative imaginary axis from T to T — i(3i. The denominator, just evolves along the 
negative imaginary axis from T io T — i(3i. The contour in the numerator may be extended 
to an arbitrary large positive time T' by inserting U{t,T')U{T' ,t) = 1 to the left of O in 



(|l^), thus becoming 

{O) (t) = Tr [U{T - zA, T)U{T, T')U{T', t)OU{t, T)] /Tr [U{T - z^, T)] (13) 

The numerator now represents the process of evolving from T < to t, inserting the operator 
O, evolving further to T', and backwards from T' to T and down the negative imaginary 
axis to T — i(3i. This process is depicted in the contour of figure (1). Eventually we take 
T — oo ; T' oo. It is straightforward to generalize to real time correlation functions of 
Heisenberg picture operators. 

This formalism allows us also to study the general case in which both the mass and the 
coupling depend on time. The insertion of the operator O, may be achieved as usual by 
introducing currents and taking variational derivatives with respect to them. 
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Because the time evolution operators have the interaction terms in them, and we would 
like to generate a perturbative expansion and Feynman diagrams, it is convenient to introduce 
source terms for all the time evolution operators in the above trace. Thus we are led to 
consider the following generating functional 



Z[J+, J-, J^] = Tr \U{T - tf3i, T; j'^)U{T, T'; J')U{T', T; .P 



(14) 



The denominator in ([T^) is simply Z[0,0,0] and may be obtained in a series expansion in 
the interaction by considering Z[0,0, J'^]. By inserting a complete set of field eigenstates 
between the time evolution operators, finally, the generating functional Z[J^ , J~ , Jl^] may 
be written as 



X 



T-ifJi 



(15) 



with the boundary conditions $+(T) = ^^{T-i(3i) = $ ; <I>+(T') = $-(T') = $2 ; $"(T) = 
$^(T) = $1. This may be recognized as a path integral along the contour in complex time 
shown in Figure (1). As usual the path integrals over the quadratic forms may be evaluated 
and one obtains the final result for the partition function 




exp 




exp 



dtC 



int \ 



^(5/5J^)|exp 



dh I dt2Jc{tl)Jc{t2)Gcih,t2] 



(16) 



Where Jc stands for the currents on the contour as shown in figure (1), Gc are the Green's 
functions on the contour and again, the spatial arguments were suppressed. 

In the limit T — 00, the contributions from the terms in which one of the currents is J"^ or 
J~ and the other is a vanish when computing correlation functions in which the external 
legs are at finite real time, as a consequence of the Riemann-Lebesgue lemma. For this real 
time correlation functions, there is no contribution from the terms that cancel between 
numerator and denominator. Finite temperature enters through the boundary conditions on 
the Green's functions (see below). For the calculation of finite real time correlation functions, 



the generating functional simplifies top8|, ^ 

Z[J+,J-] 



exp 



exp 




dt 



Cint{-iS/6J^) - Cint{iS/6J ) 



X 



T' nT' 

dti / dt2Ja{tl)Jb{t2)Gab{tl,t2) 



(17) 
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with a, 6 = +, 



The Green's functions that enter in the integrals along the contours in equations (16, 17) 
are given by (see above references) 

G++(fi,ti;f2,t2) = G>(fi,ti;f2,t2)e(ti-t2) + G'<(fi,ti;f2,t2)e(t2-ti) (18) 

G—{n,h;r2,h) = G>(fi,ti;r2,t2)e(t2-ti) + G'<(fi,ti;f2,t2)e(ti-t2) (19) 

G+-{fi,h;r2,t2) = -G<{nM\r2M) (20) 

G'-+(fi,ti;r2,t2) = -G'>(fi,ti;f2,t2) = -G<(r2,t2;ri,ti) (21) 

G'>(fi,ti;r2,t2) = {Wi.h)W2M)) (22) 

G'<(fi,T;f2,t2) = G>{n,T-il3i-r2M) (23) 

The condition (^3]) is recognized as the periodicity condition in imaginary time and is 
a result of considering an equilibrium situation for t < 0. The functions G^, are 
homogeneous solutions of the quadratic form, with appropriate boundary conditions and 
will be constructed explicitly below. 

There are many clear articles in the literature using this techniques to study real time 



correlation functions ||28||-|M 



One of our goals is to study the formation and growth of domains and the time evolution 
of the correlation functions. In particular, one relevant quantity of interest is the equal time 
correlation function 

S{r-t) = ($(f,t)<l>(0,t)) (24) 
S{r-t) = j S{k-.t) (25) 

S{k-t) = ($,^(t)$_g(t)) = -^Gt+(t;t) (26) 

where we have performed the Fourier transform in the spatial coordinates (there still is 
spatial translational and rotational invariance). Notice that at equal times, all the Green's 
functions are equal, and we may compute any of them. 

Clearly in an equilibrium situation this equal time correlation function will be time 
independent, and will only measure the static correlations. In the present case, however, 
there is a non trivial time evolution arising from the departure from equilibrium of the 
initial state. This correlation function will measure the correlations in space, and their time 
dependence. 

The function G^{ri,ti] r2,t2) will be constructed explicitly below for particular cases, in 
terms of the homogeneous solutions of the operator of quadratic fluctuations. 
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3 Rolling with fluctuations 



Inflationary cosmological models provide very appealing scenarios to describe the early evo- 
lution of our universe. Since the original model proposed by Guthfl]], several alternative 
scenarios have been proposed to overcome some of the difficulties with the original proposal. 

Among them, the new inflationary modelp[-|P is perhaps one of the most attractive. The 
essential ingredient in the new inflationary model is a scalar field that undergoes a second 
order phase transition from a high temperature symmetric phase to a low temperature broken 
symmetry phase. The expectation value (or thermal average) of the scalar field serves as 
the order parameter. Initially at high temperatures, the scalar field is assumed to be in 
thermal equilibrium and (p ^ 0. The usual field- theoretic tool to study the phase transition 
is the effective potential 0-0. At high temperatures, the global minimum of the effective 
potential is at = 0, whereas at low temperatures there are two degenerate minima. 

The behavior of the phase transition in the new inflationary model is the following: as the 
universe cools down the expectation value of the scalar fleld remains close to zero until the 
temperature becomes smaller than the critical temperature. Below the critical temperature, 
when the effective potential develops degenerate minima away from the origin, the scalar 
fleld begins to "roll down the potential hill". In the new inflationary scenario, the effective 
potential below the critical temperature is extremely flat near the maximum, and the scalar 
fleld remains near the origin, i.e. the false vacuum for a very long time and eventually rolls 
down the hill very slowly. This scenario thus receives the name of "slow rollover". During 
this stage, the energy density of the universe is dominated by the constant vacuum energy 
density Ve//(0 = 0), and the universe evolves rapidly into a de Sitter space (see for example 



the reviews by Kolb and Turner |TT[], Linde|T^ and Brandenberger [|10| ) . Perhaps the most 
remarkable consequence of the new inflationary scenario and the slow rollover transition is 



that they provide a calculational framework for the prediction of density fluct uat ions [p^j . 
The coupling constant in the typical zero temperature potentials must be flne tuned to a 
very small value to reproduce the observed limits on density fluctuations 0, |12]. 

This picture of the slow rollover scenario, is based on the static effective potential whose 
use we already criticized in the previous section. 

Guth and Pi||23|, performed a thorough analysis of the effects of quantum fluctuations on 
the time evolution. These authors analyzed the situation below the critical temperature by 
treating the potential near the origin as an inverted harmonic oscillator. They recognized 
that the instabilities associated with these upside-down oscillators lead to an exponential 
growth of the quantum fluctuations at long times and to a classical description of the proba- 
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bility distribution function. Guth and Pi also recognized that the static effective potential is 
not appropriate to describe the dynamics, that must be treated as a time dependent process. 

The purpose of this section is to study the influence of the quantum and thermal fluctu- 
ations on the "slow-rollover" picture. In particular the effect of the growth of correlations 
on the evolution of the expectation value of the infiaton field. 

In order to understand the dynamics of the phase transition and the physics of the 
instabilities mentioned above, let us consider the situation in which for time t < the 
system is in equilibrium at an initial temperature Tj > Tc where Tc is the critical temperature. 
At time t = the system is rapidly "quenched" to a final temperature below the critical 
temperature Tj < and evolves thereafter out of equilibrium. We will postpone to a later 
section the study of the evolution equations in an expanding cosmology. 

Here we just assume that such a mechanism takes place and simplify the situation by 
introducing a Hamiltonian with a time dependent mass term to describe this situation 

H{t) = |^rf=^x|itf(x) + i(V<l>(x)f + ^m2(t)<|.2(x) + A$4(^)| (27) 

m'^{t) = m^Q{-t) - iJ,^Q{t) (28) 

where both and /z^ are positive. We assume that for all times t < there is thermal 
equilibrium at temperature Tj, and the system is described by the density matrix 

p. = e-^'"' (29) 
Hi = H{t<0) (30) 

An alternative and equally valid interpretation (and the one that we like best) is that the 
initial condition being considered here is that of a system in equilibrium in the symmetric 
phase, and evolved in time with a Hamiltonian that allows for broken symmetry ground 
states, i.e. the Hamiltonian (0, ^) for t > 0. 

Using the tools of non-equilibrium statistical mechanics introduced in the previous section 
we are now in condition to obtain the evolution equations for the average of the scalar 
field in the case when the potential is suddenly changed. To account for a sudden change 
in temperature from above to below the critical temperature as described by the model 
Hamiltonian ( p7|) with (pSl) . For this purpose we use the tadpole methodQ, and write 



$±(x,t) = (/)(t) + ^±(f,t) (31) 

Where, again, the ± refer to the branches for forward and backward time propagation. The 
reason for shifting both (±) fields by the same classical configuration, is that 0(t) enters in 
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the time evolution operator as a background c-number variable, and time evolution forward 
and backwards are now considered in this background. 

The evolution equations are obtained with the tadpole method by expanding the La- 
grangian around (f){t) and considering the linear, cubic, quartic, and higher order terms in 
\E'^ as perturbations and requiring that 

< >= 0. 

It is a straightforward exercise to see that this is equivalent to one-loop to extremizing the 
one-loop effective action in which the determinant (in the log det) incorporates the boundary 
condition of equilibrium at time t < at the initial temperature. 

To one loop we find the equation of motion 



{-i)Gu{t,t) = Q 



(32) 



with Gk{t,t) = G^it,t) = G^{t,t) is the spatial Fourier transform of the equal-time Green's 
function. 

At this point, we would like to remind the reader that 



is a positive definite quantity (because the field is real) and as we argued before (and 
will be seen explicitly shortly) this Green's function grows in time because of the instabilities 



associated with the phase transition and domain growth [p^ , |23 



These Green's functions are constructed out of the homogeneous solutions to the operator 
of quadratic fluctuations 



^ + .- + M^(t) 







M\t) = {m' + ^4>^)Q{-t) + [V + ^</''W]e(t) 
The boundary conditions on the homogeneous solutions are 



uHt < 0) 

w<(/c) 



^=F«t^<{fc)t 



A 



2 



(33) 
(34) 

(35) 
(36) 



where (pi is the value of the classical field at time t < and is the initial boundary condition 
on the equation of motion. Truly speaking, starting in a fully symmetric phase will force 
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(f)i = 0, and the time evolution will maintain this value, therefore we admit a small explicit 
symmetry breaking field in the initial density matrix to allow for a small 0,. The introduction 
of this initial condition seems artificial since we are studying the situation of cooling down 
from the symmetric phase. 

However, we recognize that the phase transition from the symmetric phase occurs via 
formation of domains (in the case of a discrete symmetry) inside which the order parameter 
acquires non-zero values. The domains will have the same probability for either value of 
the field and the volume average of the field will remain zero. These domains will grow 
in time, this is the phenomenon of phase separation and spinodal decomposition familiar 
in condensed matter physics. Our evolution equations presumably will apply to the coarse 
grained average of the scalar field inside each of these domains. This average will only depend 
on time. Thus, we interpret (pi as corresponding to the coarse grained average of the field in 
each of these domains. The question of initial conditions on the scalar field is also present 
(but usually overlooked) in the slow-rollover scenarios but as we will see later, it plays a 
fundamental role in the description of the evolution. 

The identification of the initial value (pi with the average of the field in each domain is 
certainly a plausibility argument to justify an initially small asymmetry in the scalar field 
which is necesary for the further evolution of the field, and is consistent with the usual 
assumption within the slow rollover scenario. 

The boundary conditions on the mode functions U^{t) correspond to "vacuum" boundary 
conditions of positive and negative frequency modes (particles and antiparticles) for t < 0. 

Finite temperature enters through the periodicity conditions (^) and the Green's func- 
tions are 

i 1 



G>(t,t') 
G<(t,t') 



2u^{k) 1 
G>(t',t) 



(37) 
(38) 



Summarizing, the effective equations of motion to one loop that determine the time 
evolution of the scalar field are 



^+^^(t)0(t) + -03(,) + 2^(t)/^^^^coth 



■^ + .- + M^(t) 







(39) 
(40) 



with dH) , (H). 

This set of equations is too complicated to attempt an analytic solution, we will study 
this system numerically shortly. 
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However, before studying numerically these equations, one recognizes that there are sev- 
eral features of this set of equations that reveal the basic physical aspects of the dynamics 
of the scalar field. 

i) : The effective evolution equations are real. The mode functions U^{t) are complex 
conjugate of each other as may be seen from the reality of the equations, and the boundary 
conditions (|35|) . This situation must be contrasted with the expression for the effective 
potential for the analytically continued modes. 

ii) : Consider the situation in which the initial configuration of the classical field is near 



0, for t > 0, the modes for which fc^ < 



the origin (j), 
unstable. 

In particular, for early times (t > 0), when 
imately as 



ij^max) 



are 



0, these unstable modes behave approx- 



Kit) 



1 



-Wkt 



Bk = Al 



X 



fx 



Then the early time behavior of —iGk{t,t) is given by 

1 



tGk{t,t) 



1 + 



cosh{2Wkt) - 1] 



(41) 
(42) 

(43) 
(44) 



coth[Acj<(A;)/2] (45) 



This early time behavior coincides with that of the Green's function of Guth and Pi[^ and 



Weinberg and Wup2| for the inverted harmonic oscillators when our initial state (density 
matrix) is taken into account. An important point to notice is that for early times all the 
wavevectors in the unstable band grow with the same rate. 

Our evolution equations, however, permit us to go beyond the early time behavior and 
to incorporate the non-linearities that will eventually shut off the instabilities. 

These early-stage instabilities and subsequent growth of fluctuations and correlations, are 
the hallmark of the process of phase separation, and precisely the instabilities that trigger 
the phase transition. 

It is clear from the above equations of evolution, that the description in terms of inverted 
oscillators will only be valid at very early times. 

At early times, the stable modes for which k"^ > {k^axY obtained from (|41|) , ( 142|) , 



by the analytic continuation Wk — > —iu^ik) 



2 ^« 
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For t < 0, Ujl (t)Uj^ (t) = 1 and one obtains the usual result for the evolution equation 

dt^ d(f) 

with Veff{4>) the finite temperature effective potential but for t < there are no unstable 
modes. 

It becomes clear, however, that for t > there are no static solutions to the evolution 
equations for (f){t) ^ 0. 

iii) Coarsening: as the classical expectation value "rolls down" the potential hill, 
<\)^{t) increases and (kmaxit))"^ = jj? — ^(j?{t) decreases^ and only the very long-wavelength 
modes remain unstable, until for a particular time tg ; {kmax{ts)Y — 0- This occurs when 
4>'^{ts) = 2/i^/A, this is the infiexion point of the tree level potential. In statistical mechanics 



this point is known as the "classical spinodal point" and tg as the "spinodal time"|18, 19 



When the classical field reaches the spinodal point, all instabilities shut- off. From this 
point on, the dynamics is oscillatory and this period is identified with the "reheating" stage 
in cosmo logical scenarios |]12|, |10|. 

It is clear from the above equations of evolution, that the description in terms of inverted 
oscillators will only be valid at small positive times, as eventually the unstable growth will 
shut-off. 

The value of the spinodal time depends on the initial conditions of If the initial value 
(pi is very near the classical spinodal point, tg will be relatively small and there will not be 
enough time for the unstable modes to grow too much. In this case, the one-loop corrections 
for small coupling constant will remain perturbatively small. On the other hand, however, 
if (pi ^ 0, and the initial velocity is small, it will take a very long time to reach the classical 
spinodal point. In this case the unstable modes may grow dramatically making the 
one-loop corrections non-negligible even for small coupling. These initial conditions 
of small initial field and velocity are precisely the "slow rollover" conditions that are of 
interest in cosmological scenarios of "new infiation" . 

The renormalization aspects are too technical for this brief review and we refer the reader 
to our original article for details. 

As mentioned previously within the context of coarsening, when the initial value of the 
scalar field 0j ~ 0, and the initial temporal derivative is small, the scalar field slowly rolls 
down the potential hill. But during the time while the scalar field remains smaller than the 
"spinodal" value, the unstable modes grow and the one-loop contribution grows consequently. 
For a "slow rollover" condition, the field remains very small -C for a long time, 

and during this time the unstable modes grow exponentially. The stable modes, on the other 
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hand give an oscillatory contribution which is bound in time, and for weak coupling remains 
perturbatively small at all times. 

Then for a "slow rollover" situation and for weak coupling, only the unstable modes 
will yield to an important contribution to the one-loop correction. Thus, in the evolution 
equation for the scalar field, we will keep only the integral over the unstable modes in the 
one loop correction. 

Phenomenologically the coupling constant in these models is bound by the spectrum of 
density fluctuations to be within the range Xr ~ 10^^^ — 10~^^|12|. The the stable modes 



will always give a perturbative contribution, whereas the unstable modes grow exponentially 
in time thus raising the possibility of giving a non- negligible contribution. 

With the purpose of numerical analysis of the equations of motion, it proves convenient 
to introduce the following dimensionless variables 



r 



fj^Rt ; q 
Xr 



k/ fi 



R 



6yU 



^R + 1>^R<P1 



R 



(46) 
(47) 



and to account for the change from the initial temperature to the final temperature (Tj > 
Tc ; Tf < Tc) we parametrize 



fJ'R 



/^k(0) 
/ii?(0) 



- 1 



rp2 

i 
c 

1 ^ 

2^2 



(48) 
(49) 



where the subscripts (R) stand for renormalized quantities, and — /Uj:j(0) is the renormalized 
zero temperature "negative mass squared" and T"^ = 24/i|.(0)/Ai(;. Furthermore, because 
(kmaxit))'^ < and Tj > Tc, for the unstable modes Tj ^ ikmaxit)) and we can take 
the high temperature limit coth[/3jCo'<(A;)/2] 2Ti/ uo^ik). Finally the effective equations 
of evolution for t > 0, keeping in the one-loop contribution only the unstable modes as 
explained above (g^ < (gmax(T))^) become, after using uj\ = ii\{q^ + L"^) 



d^ 



dr^ 

2 



- {qmax{t)f U^{t) 
i^qmax (''") ) 

9 



+ 9'n{,r) 
= 



, ^+(r)^Y-(r) ^ 

^ 2 I r2 
q^ + L'^ 







l-3r/^(r) 
a/6X^ Ti 



27r2 



rp2 
II. 



(50) 

(51) 
(52) 
(53) 
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For Ti > Tc and Tf ^ Tc the coupling (|53| ) is bound within the range g ~ 10^^ — 
10^^. The dependence of the couphng with the temperature reflects the fact that at higher 
temperatures the fluctuations are enhanced. It is then clear, that the contribution from 
the stable modes is always perturbatively small, and only the unstable modes may introduce 
important corrections if they are allowed to grow for a long time. 

From (^) we see that the quantum corrections act as a positive dynamical renormalization 
of the "negative mass" term that drives the rolling down dynamics. It is then clear, that the 
quantum corrections tend to slow down the evolution. 

In particular, if the initial value ri{0) is very small, the unstable modes grow for a long 
time before t]{t) reaches the spinodal point i](ts) = l/y/S at which point the instabilities 
shut off. The time when this point is attained is called the spinodal time. If this is the 
case, the quantum corrections introduce substantial modifications to the classical equations 
of motion, thus becoming non- perturbative. If 7]{0) is closer to the classical spinodal point, 
the unstable modes do not have time to grow dramatically and the quantum corrections are 
perturbatively small. 

Thus we conclude that the initial conditions on the field determine whether or not the 
quantum corrections are perturbatively small. 

Although the system of equations (^, ^T]), are coupled, non-linear and integro-differential, 
they may be integrated numerically. Figures (2,3,4) depict the solutions for the classical 
(solid lines) and quantum (dashed lines) evolution (a), the quantum correction (b), i.e. the 
fourth term in ( pD]) including the coupling g, the classical (solid lines) and quantum (dashed 
lines) velocities (c) and {qmax{T))'^ (d). For the numerical integration, we have chosen 

= 1, the results are only weakly dependent on L, and taken g = 10~^, we have varied the 
initial condition r]{0) but used ^^^|t-=o = 0. 

We recall, from a previous discussion that ri{T) should be identified with the average of 
the field within a domain. We are considering the situation in which this average is very 
small, according to the usual slow-rollover hypothesis, and for which the instabilities are 
stronger. 

In figure (2. a) 77(0) = 2.3 x 10^^ we begin to see that the quantum corrections become 
important at r ^ 10 and slow down the dynamics. By the time that the classical evolution 
reaches the minimum of the classical potential at r/ = 1, the quantum evolution has just 
reached the classical spinodal i] = l/^/S. We see in figure (2.b) that the quantum correction 
becomes large enough to change the sign of the "mass term" , the field continues its evolu- 
tion towards larger values, however, because the velocity is different from zero, attaining a 
maximum (figure (2.c)) around the time when the quantum correction attains its maximum. 



17 



As 77 gets closer to the classical spinodal point, the unstability shuts off as is clearly seen in 
figure (2.d) and the quantum correction arising from the unstable modes become perturba- 
tively small. Prom the spinodal point onwards, the field evolves towards the minimum and 
begins to oscillate around it, the quantum correction will be perturbatively small, as all the 
instabilities had shut-off. Higher order corrections, will introduce a damping term as quanta 
may decay in elementary excitations of the true vacuum. 

Figures (3.a-d) show a more marked behavior for 7^(0) = 2.27 x 10"^ ; ^ = 0, notice 
that the classical evolution of the field has reached beyond the minimum of the potential at 
the time when the quantum evolution has just reached the classical spinodal point. Figure 
(3.b) shows that the quantum correction becomes larger than 1, and dramatically slows down 
the evolution, again because the velocity is different from zero (figure (3.c)) the field continues 
to grow. The velocity reaches a maximum and begins to drop. Once the field reaches the 
spinodal again the instabilities shut-off (figures (3.b,d)) and from this point the field will 
continue to evolve towards the minimum of the potential but the quantum corrections will 
be perturbatively small. 

Figures (4.a-d) show a dramatic behavior for 77(0) = 2.258 x 10^^ ; ^^^^r^ = 0. The 
unstable modes have enough time to grow so dramatically that the quantum correction 
(figure (4.b)) becomes extremely large ^ 1 (figure(4.b)), overwhelming the "negative mass" 
term near the origin. The dynamical time dependent potential, now becomes a minimum 
at the origin and the quantum evolution begins to oscillate near rj = 0. The contribution of 
the unstable modes has become non-perturbative, and certainly our one-loop approximation 
breaks down. 

As the initial value of the field gets closer to zero, the unstable modes grow for a very long 
time. At this point, we realize, however, that this picture cannot be complete. To see this 
more clearly, consider the case in which the initial state or density matrix corresponds exactly 
to the symmetric case. 77 = is necessarily, by symmetry, a fixed point of the equations of 
motion. Beginning from the symmetric state, the field will always remain at the origin and 
though there will be strong quantum and thermal fluctuations, these are symmetric and will 
sample fleld conflgurations with opposite values of the fleld with equal probability. 

In this situation, and according to the picture presented above, one would then expect 
that the unstable modes will grow indeflnetely because the scalar field does not roll down 
and will never reach the classical spinodal point thus shutting-off the instabilities. What is 
missing in this picture and the resulting equations of motion is a self-consistent treatment of 
the unstable fluctuations, that must necessarily go beyond one-loop. A more sophisticated 
and clearly non-perturbative scheme must be invoked that will incorporate coarsening, that 



18 



is the shift with time of the unstable modes towards longer wavelength and the eventual 
shutting off of the instabilities. 



4 Dynamics of Phase Separation: Domain formation 
and Growth 

In the previous section we addressed the issue of the "slow-rollover" , and we found that the 
growth of correlations that drives the process of phase separation ( "spinodal decomposition" ) 
and is the hallmark of the phase transition makes dramatic corrections to the classical 
evolution of the expectation value of the scalar field. 

Our goal in this section is to study the formation and growth of correlated domains 
and the time evolution of the correlation functions. We now consider the case of a rapid 
supercooling from the symmetric phase at Tj > Tc into the low temperature phase for Tf < Tc 
with {(f)) = 0. This value of the expectation value is a fixed point of the equations of motion. 
Thus the phase transition will occur via the formation and growth of correlated domains, 
but statistically there will be domains in which the average of the scalar field will be positive 
and others in which it will be negative, with zero average. 

In this situation, as argued in section 2 the relevant quantity of interest is the equal time 
correlation function 

S{f-t) = ($(f,t)$(0,t)) (54) 



Sir;t) = j e^'- S{k-t) (55) 

S{k-t) = («I'g(t)$_£(t)) = -^Gt+(t;t) (56) 

The function {t, t') is constructed from the homogeneous solutions to the operator of 
quadratic fluctuations 



= (57) 



with rn^it) given by (|3^). 

The boundary conditions on the homogeneous solutions are 



W±(t<0) = e^*"<(^')* (58) 

1 

uj^{k) = \k'^ + m^Y (59) 
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corresponding to positive frequency (U^(t)) (particles) and negative frequency (Uj^it)), (an 
tiparticles) respectively. 

The solutions are as follows: i): stable modes (fc^ > m'j) 



ak 
h 



— m'j 



1 / 

1 / _ uj^jkY 

2 1^ oo^ik), 



ii): unstable modes {k"^ < m?. 



U-u{t) 
W{k) 



-iaj< (k)t 



m'j — k"^ 



e(-t) + (Afce^('=)* + i?,e-^«*)e(t) 



Ak - 



Bk = {AkT 



With these mode functions, and the periodicity condition ( P5D we find 

i 1 



^^^^'^'^ 2cu<(A;)l-e-'3«-<W 
G<(t,t') = G>{t\t) 



The zeroth order equal time Green's function becomes 



2uo^{k) 



coth[/5iW<(A;)/2] 



(60) 
(61) 

(62) 
(63) 

(64) 
(65) 

(66) 
(67) 

(68) 

(69) 
(70) 

(71) 
(72) 

(73) 



for t < 0, and 



2uj^{k) 



{[1 + 2AkBk [cosh(2W^(fc)t) - 1]] e{m^f - P 



+ [1 + 2akbk [cos(2cj>(A;)t) - 1]] e{P - m^.)} coth[Acj<(A;)/2] 



(74) 
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for t > 0. 

The first term, tlie contribution of the unstable modes, reflects the growth of correlations 
because of the instabilities and will be the dominant term at long times. 



5 Zeroth order correlations: 

Before proceeding to study the correlations in higher orders in the coupling constant, it will 
prove to be very illuminating to understand the behavior of the equal time non-equilibrium 
correlation functions at tree-level. Because we are interested in the growth of correlations, 
we will study only the contributions of the unstable modes. 

The integral of the equal time correlation function over all wave vectors shows the familiar 
short distance divergences. From the above expression; however, it is clear that these may be 
removed by subtracting (and also multiplicatively renormalizing) this correlation function 
at t = 0. The contribution of the stable modes to the subtracted and multiplicatively 
renormalized correlation function is always bounded in time and thus uninteresting for the 
purpose of understanding the growth of the fluctuations. 

We are thus led to study only the contributions of the unstable modes to the subtracted 
and renormalized correlation function, this contribution is finite and unambiguous. 

For this purpose it is convenient to introduce the following dimensionless quantities 



t2 '"'i -'-i -'-c , ^ ^ /^7^-^ 

— ; L = — = — — -2 ; T = mft ; x = mjr 75 



Furthermore for the unstable modes k'^ < m'j, and for initial temperatures larger than the 
critical temperature = 24/i^/A, we can approximate coth[/3jCi;<(A;)/2] ^ 2Ti/uj^{k). Then, 
at tree-level, the contribution of the unstable modes to the subtracted structure factor (p6| 
S^^\k,t) - S^^\k,0) = (l/m/)5W(K,r) becomes 



^''\-^r) = U^U^ 1 + ^ [cosh(2iy.r)-l] (76) 




col = ^' + L' (77) 
= (78) 

To obtain a better idea of the growth of correlations, it is convenient to introduce the 
scaled correlation function 

^/ N A r'^f k^dksm(kr) , , , 

^<^-^' = ft^i l^F4r^P(fc.*)-SftO)l (79) 
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The reason for this is that the minimum of the tree level potential occurs at A$^/6mj = 1, 
and the inflexion (spinodal) point, at A$^/2mj = 1, so that V{0,t) measures the excursion 
of the fluctuations to the spinodal point and beyond as the correlations grow in time. 

At large r (large times), the product k^S{k, t) in ( ffOD has a very sharp peak ai Kg = \/t. 
In the region x < ^fr the integral may be done by the saddle point approximation and we 
obtain for Tf/T^K.^} the large time behavior 

P(x,r) ^ P(0,r)exp[-f] ^7y/f^ (80) 




Wr) - -fI^ (81) 



Restoring dimensions, and recalling that the zero temperature correlation length is ^(0) = 
we find that for T/ the amplitude of the fluctuation inside a "domain" ($^(t)), 
and the "size" of a domain ioif) grow as 

($^(t)) . ^Mf^im (S2) 

Ut) ^ (872)^^(0)^^ (83) 

An important time scale corresponds to the time at which the fluctuations of the 
field sample beyond the spinodal point. Roughly speaking when this happens, the insta- 



bilities should shut-off as the mean square root fluctuation of the field y($^(t)) is now 
probing the stable region. This will be seen explicitly below when we study the evolution 
non-perturbatively in the Hartree approximation and the fluctuations are incorporated self- 
consistently in the evolution equations. In zeroth order we estimate this time from the 
condition 3V{0,t) = 1, we use A = 10~^^ ; Tj/Tc = 2, as representative parameters (this 
value of the initial temperature does not have any particular physical meaning and was 
chosen only as representative). We find 

Ts ^ 10.15 (84) 
or in units of the zero temperature correlation length 

t ^ 14.2 ^(0) (85) 
for other values of the parameter r<j is found from the above condition on (pT]). 
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In reference |P5|| we showed explicitly that the time evolution cannot be studied in per- 
turbation theory because of the unstabilities associated with the process of phase separation 
and domain growth. These are characterized by exponentially growing contributions that 
even for small coupling will grow to become larger than the zeroth order contributions. Thus 
we need a non-perturbative treatment. 



5.1 Beyond perturbation theory: Hartree approximation 

It became clear from the analysis of the previous section that perturbation theory is in- 
adequate to describe the non-equilibrium dynamics of the phase transition, precisely because 
of the instabilities and the growth of correlations. This growth is manifest in the Green's 
functions that enter in any perturbative expansion thus invalidating any perturbative ap- 
proach. Higher order corrections will have terms that grow exponentially and faster than 
the previous term in the expansion. And even for very weakly coupled theories, the higher 
order corrections eventually become of the same order as the lower order terms. 

As the correlations and fluctuations grow, field configurations start sampling the stable 
region beyond the spinodal point. This will result in a slow down in the growth of cor- 
relations, and eventually the unstable growth will shut-off. When this happens, the state 
may be described by correlated domains with equal probabibility for both phases inside the 
domains. The expectation value of the field in this configuration will be zero, but inside each 
domain, the field will acquire a value very close to the value in equilibrium at the minimum 
of the effective potential. The size of the domain in this picture will depend on the time 
during which correlations had grown enough so that fluctuations start sampling beyond the 
spinodal point. 

Since this physical picture may not be studied within perturbation theory, we now intro- 
duce a non-perturbative method based on a self-consistent Hartree approximation. 
This approximation is implemented as follows: in the initial Lagrangian write 

A$4(r,t) = ^{<^\f,t))^\f,t) + (^$^(r-*,t) - ^(<|.^(f,t))$2(f,t)j (86) 

the first term is absorbed in a shift of the mass term 

m\t)^m'{t) + ^{^\t)) 

(where we used spatial translational invariance). The second term in (|86|) is taken as an 
interaction with the term ($^(i(:))$^(r, t) as a mass counterterm. The Hartree approximation 
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consists of requiring that the one loop correction to the two point Green's functions must be 
cancelled by the mass counterterm. This leads to the self consistent set of equations 

^'^'^^^^ = / (0 = / (0^^ i^kitPkit) coth[A^<(A:)/2] (87) 



= 



Before proceeding any further, we must address the fact that the composite operator 
($^(f, t)) needs one subtraction and multiplicative renormalization. As usual the subtraction 
is absorbed in a renormalization of the bare mass, and the multiplicative renormalization 
into a renormalization of the coupling constant. We must also point out that the Hartree 
approximation is uncontrolled in this scalar theory; it becomes equivalent to the large-N limit 
in theories in which the field is in the vector representation of 0(N) (see for example [^). 

At this stage our justification for using this approximation is based on the fact that it 
provides a non-perturbative framework to sum an infinite series of Feynman diagrams of the 
cactus type||, p^ . 

In principle one may improve on this approximation by using the Hartree propagators 
in a loop expansion. The cactus-type diagrams will still be cancelled by the counterterms 
(Hartree condition), but other diagrams with loops (for example diagrams with multiparticle 
thresholds) may be computed by using the Hartree propagators on the lines. This approach 
will have the advantage that the Hartree propagators will only be unstable for a finite time 
T < Tg. It is not presently clear to these authors, however, what if any, would be the 
expansion parameter in this case. 

It is clear that for t < there is a self-consistent solution to the Hartree equations with 
equation (p7| ) and 

($^(t)) = ($^(0-)) 

= exp[Tiuj<{k)t] (89) 



A 
2 

(90) 



colik) = k'' + m^^+^ + {^\0-)) = P + m^^j, 



where the composite operator has been absorbed in a renormalization of the initial mass, 
which is now parametrized as mjji = nWiTf /T^) — 1]. For t > we subtract the compos- 
ite operator at r = absorbing the subtraction into a renormalization of which we now 
parametrize as fn^f r = — (T|/T^^)]. We should point out that this choice of parametriza- 
tion only represents a choice of the bare parameters, which can always be chosen to satisfy 
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this condition. The logarithmic multiphcative divergence of the composite operator will be 
absorbed in a coupling constant renormalization consistent with the Hartree approximation, 
however, for the purpose of understanding the dynamics of growth of instabilities associated 
with the long-wavelength fluctuations, we will not need to specify this procedure. After this 
subtraction procedure, the Hartree equations read 

d^k 1 



mt)) - mo)) 



(27r)3 2cj<(fc) 



[U+it)U,{t)-l] coth[Aa;<(A;)/2] 



($2(0)) 



ml{t) 



T2 

- 1 

2^2 



e(-t) 



Tf 
1 ^- 

2^2 



(91) 
(92) 
(93) 



with Ti > Tc and Tf <^ Tc. With the self-consistent solution and boundary condition for 
t < 



($'(t < 0)) - ($2(0)) 
U^it < 0) 







exp[=Fia;<(/c)t] 



/c^ + mf 



iR. 



(94) 
(95) 
(96) 



This set of Hartree equations is extremely complicated to be solved exactly. However 
it has the correct physics in it. Consider the equations for t > 0, at very early times, 
when (the renormalized) (^^(t)) — ($^(0)) ~ the mode functions are the same as in the 
zeroth order approximation, and the unstable modes grow exponentially. By computing 
the expression (Rip self-consistently with these zero-order unstable modes, we see that the 
fluctuation operator begins to grow exponentially. 

As (($2(t)) — ($2(0))) grows larger, its contribution to the Hartree equation tends to 
balance the negative mass term, thus weakening the unstabilities, so that only longer wave- 
lengths can become unstable. Even for very weak coupling constants, the exponentially 
growing modes make the Hartree term in the equation of motion for the mode functions 
become large and compensate for the negative mass term. Thus when 

\r 



{^\t)) - ($2(0)) 



1 



the instabilities shut-off, this equality determines the spinodal time r^. The modes will still 
continue to grow further after this point because the time derivatives are fairly (exponen- 
tially) large, but eventually the growth will slow-down when fluctuations sample deep inside 
the stable region. 
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After the subtraction, and multiplicative renormalization (absorbed in a coupling con- 
stant renormalization), the composite operator is finite. The stable mode functions will make 
a perturbative contribution to the fluctuation which will be always bounded in time. The 
most important contribution will be that of the unstable modes. These will grow exponen- 
tially at early times and their effect will dominate the dynamics of growth and formation 
of correlated domains. The full set of Hartree equations is extremely difficult to solve, even 
numerically, so we will restrict ourselves to account only for the unstable modes. From the 
above discussion it should be clear that these are the only relevant modes for the dynam- 
ics of formation and growth of domains, whereas the stable modes, will always contribute 
perturbatively for weak coupling after renormalization. 

Introducing the dimensionless ratios ([75|) in terms of ?71/,_r ; mi,R, (all momenta are now 
expressed in units of m/^ij), dividing ( P^ by mj-j^, using the high temperature approximation 
coth[/3iCi;<(A;)/2] ^ 2Ti/uj^{k) for the unstable modes, and expressing the critical temperature 
as T'^ = 24/i/{/Aij, the set of Hartree equations ( pT] , 
differential equation for the mode functions for t > 



92) become the following integro- 



+ q^-l 



g 1^ dp- 



p" 



p"^ + L] 



■[w;(t)w-(t)-i] 



W±(t) = 



(97) 



with 



W±(t<0) = exp[Tia;<(g)t] 



9 



rp2 _ rp2 



m 



7^2 

c 



rp2 



T- 

^ 'I. 



4^2 [T2_T2]I 



(98) 
(99) 

(100) 
(101) 



The effective coupling (|101|) reflects the enhancement of quantum fluctuations by high tem- 
perature effects; for Tf/T^ ~ 0, and for couplings as weak as \r ~ 10~^^, g ~ 10~'^(Tj/Tc). 

The equations ( pT] ) may now be integrated numerically for the mode functions; once we 
flnd these, we can then compute the contribution of the unstable modes to the subtracted 
correlation function equivalent to ( [75] ) 

\r 



3P(^^)(x,r) 



(<l>(f,t)<l>(0,t))-(<l>(f,0)$(0,0)) 
sin(px) 



p- 



px 



u;{t)u;it)-i 



(102) 
(103) 
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In figure (4) we show 

-^(($2(r)) - ($2(0))) = 3(P^^(0, r) - V^^{0, 0)) 

(solid fine) and also for comparison, its zeroth-order counterpart 3(I?(°)(0,r) - D(°)(0,0)) 
(dashed line) for Aj? = 10~^^ , Ti/Tc = 2. (Again, this value of the initial temperature does 
not have any particular physical significance and was chosen as a representative). We clearly 
see what we expected; whereas the zeroth order correlation grows indefinitely, the Hartree 
correlation function is bounded in time and oscillatory. At r 10.52 , 3{'D^^^\Q,t) — 
T>'^^^\0,t)) = 1, fiuctuations are sampling field configurations near the classical spinodal, 
fiuctuations continue to grow because the derivatives are still fairly large. However, after this 
time, the modes begin to probe the stable region in which there is no exponential growth. 
At this point ^j^f (($^(r)) — $^(0))), becomes small again because of the small couphng 

f,R 

g 10~^, and the correction term diminishes. When it becomes smaller than one, the 
instabilities set in again, modes begin to grow and the process repeats. This gives rise to an 
oscillatory behavior around (('l'^(r)) — $^(0)) = 1 as shown in figure (5). 

f,R 

In figures (6.a-d), we show the structure factors as a function of x for r = 6, 8, 10, 12, 
both for zero-order (tree level) T>^^^ (dashed fines) and Hartree (solid fines). These 

correlation functions clearly show that correlations grow in amplitude and that the size of 
the region in which the fields are correlated increases with time. Clearly this region may 
be interpreted as a "domain" , inside which the fields have strong correlations, and outside 
which the fields are uncorrelated. 

We see that up to the spinodal time 10.52 at which 2^((^'^ ('^^)) ~ ^^(0))) = 1' 
the zeroth order correlation SV^^^Ojt) is very close to the Hartree result. In fact at r^, 
the difference is less than 15%. In particular for these values of the coupling and initial 
temperature, the zeroth order correlation function leads to Tg ~ 10.15, and we may use the 
zeroth order correlations to provide an analytic estimate for r^, as well as the form of the 
correlation functions and the size of the domains. 

The fact that the zeroth-order correlation remains very close to the Hartree-corrected 
correlations up to times comparable to the spinodal is clearly a consequence of the very 
small coupling. 

To illustrate this fact, we show in figures (7,8) the same correlation functions but for 
A = 0.01, Ti/Tc = 2. Clearly the stronger coupling makes the growth of domains much 
faster and the departure from tree-level correlations more dramatic. Thus, it becomes clear 
that for strong couplings, domains will form very rapidly and only grow to sizes of the order 
of the zero temperature correlation length. The phase transition will occur very rapidly, and 
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clearly our initial assumption of a rapid supercooling will be unjustified. This situation for 
strong couplings, of domains forming very rapidly to sizes of the order of the zero temperature 
correlation length, is the picture presented by Mazenko and collaborators However, for 
very weak couplings (consistent with the bounds from density fluctuations), our results 
indicate that the phase transition will proceed very slowly, domains will grow for a long time 
and become fairly large, with a typical size several times the zero temperature correlation 
length. In a sense, this is a self consistent check of our initial assumptions on a rapid 
supercooling in the case of weak couplings. 

Thus, as we argued above, for very weak coupling, we may use the tree level result to give 
an approximate bound to the correlation functions up to times close to the spinodal time 
using the result given by equation (pT|), for Tj ^ 0. Thus, we conclude that for large times, 
and very weakly coupled theories {Xr < 10~^^) and for initial temperatures of the order of 
the critical temperature, the size of the domains ^oit) will grow typically in time as 



with ,^(0) the zero temperature correlation length. The maximum size of a domain is 
approximately determined by the time at which fluctuations begin probing the stable region, 
this is the spinodal time Tg and the maximum size of the domains is approximately C,d{'^s)- 

An estimate for the spinodal time, is obtained from equation (|8TD by the condition 
3V[ts) = 1, then for weakly coupled theories and Tj ^ 0, we obtain 



In 



3A V i§-f 



/ ^2 i 



(105) 



It is remarkable that the domain size scales as C,D{t) ~ just like in classical theories 
of spinodal decomposition, when the order parameter is not conserved, as is the case in 
the scalar relativistic field theory under consideration, but certainly for completely different 
reasons p7[|. At the tree level, we can identify this scaling behavior as arising from the 
relativistic dispersion relation, and second order time derivatives in the equations of motion. 



a situation very different from the classical description of the Allen-Cahn-Lifshitz[^, IIP 
theory of spinodal decomposition based on a time- dependent Landau- Ginzburg model. 

When the scalar coupling is strong, the phase transition proceeds rapidly, and domains 
will not have time to grow substantially. Their sizes will be of the order of the zero temper- 
ature correlation length. This will become important in a later section when we study the 
dynamics in a phenomenologically relevant strongly coupled theory. 
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Beyond Hartree: 

The Hartree approximation, keeping only the unstable modes in the self- consistent equa- 
tion, clearly cannot be accurate for times beyond the spinodal time. When the oscillations in 
the Hartree solution begin, the field fluctuations are probing the stable region. This should 
correspond to the onset of the "reheating" epoch, in which dissipative effects become im- 
portant for processes of particle and entropy production. Clearly the Hartree approximation 
ignores all dissipative processes, as may be understood from the fact that this approximation 
sums the cactus type diagrams for which there are no multiparticle thresholds. Furthermore, 
in this region, the contribution of the stable modes to the Hartree equation becomes impor- 
tant for the subsequent evolution beyond the spinodal point and clearly will contribute to 
the "reheating" process. A possible approach to incorporate the contribution of the stable 
modes may be that explored by Avan and de Vega[^| in terms of the effective action for the 
composite operator. 

Thus, although the Hartree approximation may give a fairly accurate picture of the 
process of domain formation and growth, one must go beyond this approximation at times 
later than the spinodal time, to incorporate dissipative effects and to study the "reheating" 
period. Clearly, one must also attempt to study the possibility of "percolation of domains" . 
Furthermore, the Hartree approximation is essentially a Gaussian approximation, as the 
wave-functional (or in this case the functional density matrix) is Gaussian with kernels that 
are obtained self- consistently. The wave-functional must include non-gaussian correlations 
that will account for the corrections to the Hartree approximation and will be important to 
obtain the long time behavior for t > Tg. 

6 Non-Equilibrium Dynamics in FRW Cosmologies 

In this section we provide an alternative and equally powerful method to obtain the evo- 
lution equations by evolving in time an initially prepared density matrix. This method is 
fundamentally equivalent to the one described in the previous section, but is restricted to 
gaussian density matrices. This restriction, however still offers the possibility of studying 
both perturbative and non-perturbative approximations and we include it here as one more 
tool in non-equilibrium field theory. Furthermore this method permits the study of effects 
of particle production via parametric amplification in a more natural way. 

What we address in this section is the effective evolution equations for the order parameter 
and quantum fluctuations in homogeneous and isotropic FRW expanding cosmologies. 

By and large, the various models of inflation make the assumption that the dynamics 
of the spatial zero mode of the (so-called) inflaton field is governed by some approximation 
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to the effective potential which incorporates the effects of quantum fluctuations of the field. 
Thus the equation of motion is usually of the form: 



(^ + 3-0 + Kff(0)=O 
a 



(106) 



The problem here is that the effective potential is really only suited for analyzing static 
situations; it is the effective action evaluated for a field configuration that is constant in 
time. Thus, it is inconsistent to use the effective potential in a dynamical situation. Notice 
that such inconsistency appears for any inflationary scenario (old, new, chaotic, ...). 

Our approach is to use the functional Schrodinger formulation, wherein we specify the 
initial wavefunctional \I^[$(7);t] (or more generally a density matrix p[$(T^, $(7); t]), and 
then use the Schrodinger equation to evolve this state in time. We can then use this state to 
compute all of the expectation values required in the construction of the effective equations 
of motion for the order parameter of the theory, as well as that for the fluctuations. The 
Schrodinger approach has already been used in the literature at zero temperature [^] and to 



study non-equihbrium aspects of field theories|41, 42 



We start by setting up the Schrodinger formalism for spatially flat FRW cosmologies. 
Consider a scalar field in such a cosmology where the metric is: 



ds"^ = dt^ - a^{t)dx^ 
The action and Lagrangian density are given by 
S = I d'^xC 



a^it) 
1 



1a2/- n l(V$(x,t))2 

2 ^ ' ^ 2 a(t)2 V V > 

A 



7^ 



6 - + ^ 
V a 



with 7^ the Ricci scalar. The canonical momentum conjugate to $ is 

n(x,t) = a^{t)^{x,t) 



and the Hamiltonian becomes 
E(t)- 



d^x { — r-- 



a(t) 



(V<l>)2 + a3(t)F($) 



(107) 

(108) 
(109) 

(110) 
(111) 

(112) 
(113) 
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In the Schrodinger representation (at an arbitrary fixed time to), the canonical momentum 
is represented as 

n(f) = -ih- 



Wave functional obey the time dependent functional Schrodinger equation 

in^^^^^H^ij[<^,t] (114) 

Since we shall eventually consider a "thermal ensemble" it is convenient to work with a 
functional density matrix p with matrix elements in the Schrodinger representation p[$(7), $(7); i]. 
We will assume that the density matrix obeys the functional Liouville equation 

th^ = [H{t),p] (115) 

whose formal solution is 

p{t)^U{t,to)p{to)U-\t,to) 

where U{t,to) is the time evolution operator, and p{to) the density matrix at the arbitrary 
initial time to- 

The diagonal density matrix elements t] are interpreted as a probability density 

in functional space. Since we are considering an homogeneous and isotropic background, the 
functional density matrix may be assumed to be translationally invariant. Normalizing the 
density matrix such that Trp — 1, the order parameter is defined as 

Ht) = ^/ d'xmx,t)) = d^xTr [pmi^)] = ^/ d'xTr [p{to)U-\t,toMx,to)U{t,to) 

(116) 

where Q is the comoving volume, and the scale factors cancel between the numerator (in the 
integral) and the denominator. Note that we have used the fact that the field operator does 
not evolve in time in this picture. The evolution equations for the order parameter are the 
following 

^ - ^/^'^<n,^,0> ^ ^j'e^Trmni.l ^ ^ (UT) 

'-^ = / (118) 

It is now convenient to write the field in the Schrodinger picture as 

^[x) = (f){t) +ri{x,t) (119) 
{7j{x,t)) = (120) 
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Expanding the right hand side of ( |118| ) we find the effective equation of motion for the 
order parameter: 

^'^W+3^^4^ + nW) + n^ /*':^{.)'(f.0) + - = (121) 



df^ a{t) dt '"^^ 2n 

where primes stand for derivatives with respect to 0. To leading order in the loop expansion 
we need that ijf{x^t)) = 0{h). This will be guaranteed to this order if the density matrix 
is assumed to be Gaussian with a covariance (width) 0{l/h). If ( |119|) is introduced in the 
Hamiltonian, we arrive at: 

m = J d'x {-^^ + ^ (Vr/)' + a%t) (Vi<p) + VWv + \v"{c^W + 

(122) 

Keeping only the terms quadratic in rj in (|122|) gives the first order term in the loop 
expansion. 

It is convenient to introduce the discrete Fourier transform of the fields in the comoving 
frame as 

Vix,t) = -^Y.vdt)e-''' (123) 

In this representation, the quadratic approximation to the Hamiltonian ( |122D becomes the 
Hamiltonian of a collection of independent harmonic oscillators for each mode k 

oolit) = a{t)P + a''{t)V"{<P{t)) (125) 

We propose the following Gaussian ansatz for the functional density matrix elements in 
the Schrodinger representation 



+^vrfc(t)(r7_,(t)-77_fc(t))| (126) 

Vk{t) = $fc - 0(t)v^5g^o (127) 
Ut) = $fc-0(t)v^5fco (128) 
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where = and TCkit) is the Fourier transform of (n(x)). This form of the density 

matrix is dictated by the hermiticity condition p'^[$,$,t] = as a result of this, 

Bk{t) is reaL The kernel Bk{t) determines the amount of "mixing" in the density matrix, 
since if -Bfc = 0, the density matrix corresponds to a pure state because it is a wave functional 
times its complex conjugate. 

In order to solve for the time evolution of the density matrix ( |115|) we need to specify 
the density matrix at some initial time to- It is at this point that we have to assume 
some physically motivated initial condition. We believe that this is a subtle point that 
has not received proper consideration in the literature. A system in thermal equilibrium 
has time-independent ensemble averages (as the evolution Hamiltonian commutes with the 
density matrix) and there is no memory of any initial state. However, in a time dependent 
background, the density matrix will evolve in time, departing from the equilibrium state and 
correlation functions or expectation values may depend on details of the initial state. 

We will assume that at early times the initial density matrix is thermal for the modes 
that diagonalize the Hamiltonian at (we call these the adiabatic modes). The effective 
temperature for these modes is ksTo = l/f3o- It is only in this initial state that the notion of 
"temperature" is meaningful. As the system departs from equilibrium one cannot define a 
thermodynamic temperature. Thus in this case the "temperature" refers to the temperature 
defined in the initial state. 

The initial values of the order parameter and average canonical momentum are 0(to) = 0o 
and nito) = tTo respectively. Defining the adiabatic frequencies as: 



+ V"{(t)cl{to)) 



(129) 



we find that the initial values of the time dependent parameters in the density matrix (|126| ) 



are 



Ak{to) 
Bk{to) 

4>{to) 



AUQ = WkiQa'iQcothiPohWkit,)] 
sinh [PonWkiQ] 



7rh 

vr(t. 



tanh 



Tin 



(130) 
(131) 

(132) 
(133) 



The initial density matrix is normalized such that Trp{to) = 1. Since time evolution is 
unitary such a normalization will be constant in time. For To = the density matrix 
describes a pure state since Bk = 0. 
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As an example, consider the case of de Sitter space. The scale factor is given by a(t) = 
ttoC^^ and for To 0, to — > — oo we recognize the ground state wave-functional for the 
Bunch-Davies vacuum[^, For To 7^ this initial density matrix corresponds to a 

thermal ensemble of Bunch-Davies modes. Certainly this choice is somewhat arbitrary but 
it physically describes the situation in which at very early times the adiabatic modes are in 
local thermodynamic equilibrium. Whether or not this situation actually obtains for a given 
system has to be checked explicitly. In the cosmological setting, the nature of the initial 
condition will necesarily have to result from a deeper understanding of the relationship 
between particle physics, gravitation and statistical mechanics at very large energy scales. 

Although we will continue henceforth to use this thermal initial state, it should be em- 
phasized that our formalism is quite general and can be applied to any initial state. 

In the Schrodinger picture, the Liouville equation (|115|) becomes 



ih 



dp[<^,^,t] 
dt 



E 



6^ 



6^ 



2a^{t) \6rik6r]_k Srik6r]_k, 



W{t)V'^k{(p{t)) iVk - Vk) + 2^fc(^) (VkV^k - VkV-k 



(134) 



Since the modes do not mix in this approximation to the Hamiltonian, the equations for the 
kernels in the density matrix are obtained by comparing the powers of rj on both sides of the 
above equation. We obtain the following equations for the coefficients: 



.A 

iAk 
iBk 



1 



2a3(t) 



ait) 



(Ak - Al) 



-{Ak-Al 



am 



TC 



a^{t) 



(135) 
(136) 

(137) 

(138) 
(139) 



The last two equations are identified with the classical equations of motion for the order 
parameter ( p,21| ). The equation for Bk{t) refiects the fact that a pure state Bk = remains 
pure under time evolution. 

Writing Ak in terms of its real and imaginary components Ak(t) = A^kit) + iAik(t) (and 
because Bk is real) we find that 

Bkit) Bk{to) 



ARkit) Ankito) 



(140) 
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and that the time evolution is unitary (as it should be), that is 

A4(t) 



constant 



(141) 



The initial conditions ( |130| , |131| ) and the invariance of the ratio (|14CI|) suggest that the 
solution for the real part of A and for B may be obtained by introducing a complex function 
Ak{t) = ABkit) + iAik{t) 



ARk{t) = ARk{t) tanh IPohWkito)] = -Bk{t) sinh [(3ohWk{to)\ 
ArM = WM a\Q 
Aik{t) = Ajk{t) 

In this form, the real and imaginary parts of A satisfy the equations 

2 



ARk{t) 

-Aik{t) 



a3(t) 
1 



ARkit)Aikit) 



(142) 
(143) 
(144) 



(145) 
(146) 



These two equations may be combined in one complex equation for the combination Ak{t) 
ARk(t) + iAik{t) that obeys the Riccati-type equation 

1 



iAkit) 



with the initial conditions: 



ARkito) 

Aik{to) 



Ak'{t)-ujt{t)a%t) 



Wk{to)a\to) 




(147) 



(148) 
(149) 



The Riccati equation ( |147| ) becomes a more amenable differential equation by the change 
of variables 



Akit) = -ia\t) 



(150) 



The solution to the Ricatti equation with the above initial conditions is detailed in 
appendix A. We find that it is convenient to introduce two real mode functions (for each 
wavevector k) and write 



<^k{t) 



Ukljt) + ^k2{t) 

'a\t)Wk{to) ' 



(151) 
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These mode functions obey a Schrodinger-like equation. The initial conditions on ipk(to) are 
(see appendix A) 



'-Pkito) 



■^Rkit) is given by equation ( |215| ) in appendix A, so that: 



ARk{t) 

Bk{t) 



1 



1 



COth [(^JlWkito) 

1 



The equal time two-point function for the fluctuation becomes 

n h 



{Vk{t)v-k{t)) 



Mt) rcoth[l3ohWkito)/2] 



2{ARk{t) + Bk{t)) 2' 
The one-loop equation of motion for the order parameter thus becomes 



(152) 
(153) 

(154) 



^ + 3-<P + V\<P) + V'"{<f>)^ 



coth[pohWk{to)/2] = 



(155) 



(27r)3 2 

with the function (pk{t) defined in appendix A by (|208| , |209|) in which, to this order in h, 
only the classical solution (f)d{t) enters. A consistent numerical solution of these equations 
to 0{h) would involve splitting (j){t) = (f)ci{t) + ^4'{i){t) and keeping only the 0{h) terms in 
the evolution equation. This will result in two simultaneous equations, one for the classical 
evolution of the order parameter and another for 0(i)(t). 

This equation of motion is clearly very different from the one obtained by using the 
effective potential. It may be easily seen (by writing the effective action as the classical 
action plus the logarithm of the determinant of the quadratic fluctuation operator) that this 
is the equation of motion obtained by the variation of the one-loop effective action. 

The static effective potential is clearly not the appropriate quantity to use to describe 
scalar field dynamics in an expanding universe. Although there may be some time regime 
in which the time evolution is slow and fluctuations rather small, this will certainly not be 
the case at the onset of a phase transition. As the phase transition takes place, fluctuations 



become dominant and grow in time signaling the onset of long range correlations |T6|, |35 . 
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6.1 Hartree equations 

Motivated by our studies in Minkowski space reported in the previous section, in which 
showed that the growth of correlation and enhancement of fluctuations during a phase tran- 
sition may not be described perturbatively, we now proceed to obtaining the equations of 
motion in a Hartree approximation. As argued in the previous section, this approximation 
is non-perturbative in the sense that it sums up infinitely many diagrams of the cactus- 
type@, ^. 



Vilenkin||44]| has previously studied a simplified version of the Hartree approximation. 
The Hartree self-consistent approximation is implemented as follows. We decompose the 
field as in ( |119D , using a potential as in ( |11(J| ). We find that the Hamiltonian becomes 



+ ^W + ^Ar/^H (156) 



3! ^ ' 4! 

The Hartree approximation is obtained by assuming the factorization 

r]^{x,t) ?>{rf{x,t))r]{x,t) (157) 
r/^(f,t) Q{^f{x,t))^f{x,t)-^rf{x,t))'^ (158) 

where (■ ■ ■) is the average using the time evolved density matrix. This average will be deter- 
mined self-consistently (see below). Translational invariance shows that {jf{x,t)) can only 
be a function of time. This approximation makes the Hamiltonian quadratic at the expense 
of a self-consistent condition. In the time independent (Minkowski) case this approxima- 
tion sums up all the "daisy" (or "cactus") diagrams and leads to the self-consistent gap 
equationp|. In this approximation the Hamiltonian becomes 

{"2^1^ + ^ {^vf + a%t) (V«(0), + lv^'\<^W) } (159) 
V(0) = V(0) - hivr (160) 
V«(0) = ^'(0) + ^A0(r^^) (161) 
V^^\<P) = V"i<P) + h{r^') (162) 
We can now introduce the Fourier transform of the field as in (11 231). The Hamiltonian will 



have the same form as ( |124| ) but the time dependent frequencies ( |125| ) and linear term in rj 
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become 



uKt) = a{t)k^ + a%t)V^'\(P{t)) (163) 
(0(t)) = V^'Hm^h. (164) 



The ansatz for the Gaussian density matrix is the same as before (|126|) , as are the evolution 
equations for the coefficients A^it), Bk{t), Afkit). However, the frequencies are now given 
by ( |163| ). The c/asszca/ equations of motion ( p.38|jl39| ) now become 

_vr = V^^\<P{t))a^{t) (165) 

(166) 



a3(t) 

The equations for the coefficients Afc(t), Bk{t), N'kit) are again solved in terms of the mode 
functions given in appendix A but with V""(0c/(i)) now replaced by V^'^\(p{t)) and the fol- 
lowing replacement of the adiabatic frequencies: 



a{to) 



(167) 
(168) 



The mode functions of the appendix obey the differential equations and initial conditions 
with this replacement and 

I Vk{t) (t) P . 

The initial conditions at to are now given in terms of these new adiabatic frequencies (see 
appendix A). The equal time two-point function thus becomes 

h r d^k 



(169) 



which leads to the following set of self-consistent time dependent Hartree equations: 



d'k I (t) 
(27r)3 2 



coth[(3ohWkito)/2] = 



h f d'k I (t) 

(27r)3 2 



(170) 
(171) 
(172) 
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m{to) 

\ a' (to) 



(173) 



(174) 



The renormalization aspects of these equations have been worked out in detail in reference [p£ 
to which we refer the reader. 



We find||3^ that the renormahzed self-consistent Hartree equations become after letting 
A = oo: 



. + 3^0 + + Cnn4> + + ^{v')r = 
d'^ a{t) d k'^ 2 , c ^ , a2 , -^R / 2\ 



(175) 



R 



vUt) = (176) 



where (^7^)_r is the renormalized composite operator 

For completeness we quote the renormalized equation of motion for the order parameter 
to one-loop 



R 



rf2 a(t) d 



(177) 

(178) 
(179) 



The one-loop coupling constant renormalization differs from that in the Hartree approxi- 
mation by a factor of three. This is a consequence of the fact that the Hartree approximation 
is equivalent to a large N approximation and only sums up the s-channel bubbles. 



6.2 High Temperature Limit 

One of the payoffs of understanding the renormalization aspects and the large momentum 
behaviour of the various fluctuation integrals is that that it permits the evaluation of the 
high temperature limit. We shall perform our analysis of the high temperature expansion 
for the Hartree approximation. The one-loop case may be read off from these results. 
The finite temperature contribution is determined by the integral 



J = h 



d'k I (t) 



(180) 
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For large temperature, only momenta k > To contribute. The leading contribution is found 

1 



to bef39 



J 



I2h 



kBToa{to) 



a{t) 



[1 + Oil/To) + ■ 



181) 



Thus we see that the leading high temperature behavior reflects the physical red-shift in 
the cosmological background and it results in an effective time dependent temperature 



Teffit) = To 



a{t) 



To leading order, the expression obtained for the time dependent effective temperature 
corresponds to what would be obtained for an adiabatic (isentropic) expansion for blackbody- 
type radiation consisting of massless relativistic particles evolving in the cosmological back- 
ground. 

This behavior only appears at leading order in the high temperature expansion. There 
are subleading terms that must be taken into account. These can be calculated within the 
high temperature expansion and we do this below. To avoid cluttering of notation, we will 
set ks = h = 1 in what follows. We find the linear and logarithmic contributions in To to 
be|| 



Ji 



llin 



a(to) 
a{t) 



Tom{To) 



dz{a^(t)m(To)\fz^ 



Note that the above integral is finite. 

We obtain after some straightforward algebra PU[: 



a{t) 
2a{to) 



;i82) 



Ji 



In 



llog 



m{To) 
To 



n 

' 6" 



+ 



a{to) 


2 - 


a{t) 





m\To) + 



MI 

24 



MS 

24 



;i83) 



That is, in the limit To » m{To), Ji = Jmn + Jiiog + 0((To)°). 

Comparing the 0{T^, To, In To) contributions it becomes clear that they have very differ- 
ent time dependences through the scale factor a{t). Thus the high temperature expansion as 
presented will not remain accurate at large times since the term quadratic in To may become 
of the same order or smaller than the linear or logarithmic terms. The high temperature 
expansion and the long time limit are thus not interchangeable, and any high temperature 
expansion is thus bound to be valid only within some time regime that depends on the initial 
value of the temperature and the initial conditions. 
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As an illustration of this observation, we calculate Jmn explicitly in the case of de Sitter 
space. We need to obtain | (p^(t) p in order to evaluate the integral in (|182|) . Inserting 
the term proportional to in the Hartree equations, we find that | ^^i^) P obeys the 
differential equation 



24 



MS 

24 



The solution of this equation is given by: 



;i84) 



Bk 



r(2), 



1 

H 



24 



(185) 



fl 2") 

where -^3/2 the Hankel functions and we have assumed 



24 



The coefficients Ci, C2 are determined by the initial conditions on ^p^i^) described above. 
We finally obtain 



llin 



1 z 



^ + w?{T)z 



24 



(186) 



This term is time independent, finite and positive. This example clearly illustrates the fact 
that different powers of To enter in the expansion with different functions of time and that 
the high temperature expansion is non-uniform as a function of time. 



6.3 Evolution of the Initial Distribution 



The initial density matrix at t = to was assumed to be thermal for the adiabatic modes. 
This corresponds to a Boltzmann distribution for the uncoupled harmonic oscillators for the 
adiabatic modes of momentum fc, and frequencies VVfc(to) (in the Hartree approximation; as 
usual, the one-loop result can be found by replacing this with H4(to))- That is 



p{to) 



(187) 
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The creation and annihilation operators define the initial occupation number of the adiabatic 
modes: 



"I(^o)afc(tc 
1 



(189) 
(190) 



where the expectation value in (|19CI|) is in the initial density matrix at time to- 

In a time dependent gravitational background the concept of particle is ill-defined. How- 
ever, by postulating an equilibrium initial density matrix of the above form, a preferred 
"pointer" basis is singled out at the initial time. It is this basis that provides a natural 
definition of particles at the initial time and we can use it to ask: how does the expectation 
value of this number operator evolve in time? 

At any time t, this expectation value is given by 



Tr 



a\{to)ak{to)p{t) 



Trpito) 



(191) 



This quantity gives information on how the original Boltzmann distribution function for the 
adiabatic modes evolves with time. The /c = mode will receive a contribution from the 
order parameter, but since the number of particles is not conserved (no charge) there is no 
bose condensation and the = mode will give a negligible contribution to the total number 
of particles. Thus we only concentrate on the A; 7^ modes. 

The expectation value ( [191| ) may be easily computed by writing the creation and anni- 
hilation operators in terms of 77^; = 6/6r]_k in the Schrodinger picture at to- The result 
of doing this is: 



2h 
1 



U ^ 1 

al{to) = 
akito) = 

After some straightforward algebra we find 



+ a\to)Wk{to)v~k 
+ ^a^{to)m{to)Vk 



{Nk){t) + - 

I ^kit, to) 1^ 



2 I Mt,to) P -1) i{Nk)ito) + - 



1 + 



1 ^^{t) 


2 " 


<(0) 


2 



a%t) I 0f (t) 
a%to)Wi{to) I <(t) 



+ 



(192) 
(193) 



(194) 
(195) 
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where we have made use of ( pi5| , pi6| ). This result exhibits the two contributions from 
"spontaneous" (proportional to the initial thermal occupation) and "induced" (independent 
of it). 

We now show that this result may be understood as a Bogoliubov transformation. To do 
this, consider the expansion of the field in the Heisenberg picture: 



(196) 



where the mode functions satisfy 



^fe + 3 — ^ + 



a dt 



+ v(^)(0W) 



(197) 



together with the self-consistency relation. Then the Heisenberg field ?7fc(t) is a solution of 
the Hartree Heisenberg equations of motion and the al , ak create and destroy the Hartree- 
Fock states. Notice that in the Heisenberg picture, these creation and annihilation operators 
do not depend on time. Using the Wronskian properties of the functions ipk{t) (see appendix 
A) we can invert and find the creation and annihilation operators in terms of T]k{t) and its 
canonically conjugate momentum H_fc(t). Once we have expressed these operators in the 
Heisenberg picture in terms of the field and its canonically conjugate momentum, we can go 
to the Schrodinger picture at time to- In this picture the creation and annihilation operators 
depend on time and are given by 



ak{t) 
al,(t) 



—i r 
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nk{to)^'^\t) - a\t)r^u{to)^u\t) 



(198) 
(199) 



The Schrodinger picture fields at to can be written in terms of the operators 



and we finally find the creation and destruction operators at time t to be related to those at 
time to by a Bogoliubov transformation: 



ttfc(t) = T+^k{t,to)ak{to) + J^-,k(t,to)a'ik(to) 



(200) 



If we now compute the average of the new creation and annihilation operators in the initial 
density matrix and write the mode functions ipk{t) in terms of the real functions Ui^2 as 
defined in appendix A, we recognize | J^+^k(t,to) P to be the same as | J-'k(t,to) P given by 

(HI). 

We also find that 



^ +,k(t, to) 



^-,k{t, to) \^— 1 
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as is required for a Bogoliubov transformation. 

One way to interpret this result is that, at least within the one-loop or Hartree approxi- 
mations, time evolution corresponds to a Bogoliubov transformation. This interpretation is, 
in fact, consistent with the result that in these approximation schemes, the density matrix 
remains Gaussian with the only change being that the covariance and mixing terms change 
with time. 

Thus within the one-loop or Hartree approximation, time evolution corresponds to a 
"squeezing" of the initial state. The covariance changes with time and this corresponds to 
a Bogoliubov transformation. As argued by Grishchuk and Sidorov[^] the amplification of 
quantum fluctuations during inflation is a process of quantum squeezing and it corresponds 
to a Bogoliubov transformation. The properties of these "squeezed" quantum states have 
been investigated in references ^ 



Hu and Pavon||49[], Hu and Kandrup0] and Kandrup|5T|] have introduced a non-equilibrium 



coarse-grained entropy that grows in time as a consequence of particle production and "para- 
metric amplification". This definition was generalized by Brandenberger et al. [|^, and 



Gasperini and collaborators to give a measure of the entropy of the gravitational field. 
The growth of this entropy is again a consequence of the parametric amplification of fluctu- 
ations and the "squeezing" of the quantum state under time evolution. 

These authors argue that the non-equilibrium coarse-grained entropy in the mode of 
(comoving) wavevector k is 

Sk^\n{{Nk)it)) (201) 

in the case when {Nk){t) ^ 1. Thus the growth of entropy is directly associated with 
"particle production" or in our case to the evolution of the initial Boltzmann distribution 
function. 

The coefficient of parametric amplification is related to the Bogoliubov coefficient given 
by equation ( p.95| ). Thus this coefficient directly determines the time dependence of the 
non-equilibrium coarse-grained entropy. 

6.4 Appendix A 



The Riccati equation (|147|) can be transformed into a linear differential equation by the 
change of variables 



= -ta'it)^ (202) 
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We find that ipk(t) obeys a simple evolution equation 



a dt 



(203) 



with (pciit) the solution to the classical equation of motion ( |138| , |139| ). In the Hartree ap- 
proximation V"{4)ci{t)) should be replaced by V'-^^(0(t)) with the full solution of the 
self-consistent equations. 



dt^ 



a dt 



(204) 



The Wronskian for two arbitrary solutions to the above differential equation is 

C 



a3(t) 



(205) 



with C a constant. Then writing ifk{t) = fik{t) + i^2k{i) with ipi 2 real solutions we find 

C 



ARk{t) 

Aik{t) 



¥>lk + Vlk 

3/^^ ^'ikVlk + V2kV2k 



-a\t) 



Vlk + ^Ik 



(206) 
(207) 



It proves convenient to introduce (for all k) the real functions Wi,2(t) as 

<^l,2(t) = [a{t)] 2. 



Wk{to) 



(20^ 



The Uak for a = 1, 2 are real and satisfy the Schrodinger-like differential equation 

d^ 



3 f d 1 k"^ T./// , N / N 



Kak{t) = 



W[U^^kM2,k] = CWkit, 



(209) 
(210) 



with W[- ■ •] the Wronskian, and C is the same constant as above. Since the choice of C 
corresponds to a choice of normalization of these functions, we choose C = 1. The initial 
conditions ( p.48| , |149|) still leave one free condition on these functions, we choose 



Uikito) ^0 

U2k{to) = 



(211) 
(212) 
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The boundary conditions on the mode functions Ua,k are 

3a(to; 



2a{to) 



U2k{to) = Wu{to) 



a? [to] 



+ V"{ct)a{U)) 



(213) 
(214) 



and the corresponding replacement for the Hartree case. Thus the final solution to the 
Riccati equation (|147|) with the given initial conditions is 



ARkit) 

Aikit) 



-a\t) 



Wife (^Ifc — it^lfc) + ^2A: (^2fc — if ^2fc 



(215) 
(216) 



In terms of the original functions ^k{t) ( P02|) the initial conditions are simply 

^k{to) = , ^ (217) 



a%QWk{to) 



(218) 



Then the initial conditions for ipk{t) ', fli't) are naturally interpreted as those for negative 
and positive (adiabatic) frequency modes at the initial time to- 

The reason for introducing the functions Ua,k{t) is because these obey a simpler second 
order Schrodinger-like equation which is amenable to be studied in the asymptotic regime 
via WKB approximations (see the section on renormalization). 



6.5 Appendix B: Conformal Time Analysis 

It is interesting to see how some of our results can be obtained by rewriting the metric in 
terms of the conformal time defined by: 

The first thing we should note is that the physics should not depend on what time 
coordinate is used, since the theory should be generally coordinate invariant (there are no 
gravitational anomalies in four dimensions). Thus, the field amplitude and the canonical 
momentum in conformal time should be related to those in comoving time via a canonical 
transformation. We now show that this is indeed the case. 
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Using the conformal time version of the hne element ds^ = a?{ri){dif' — dx^), the scalar 
field action can be written as: 



= J dr]d^xa\r]) 



1 



2a2(r7) 



((9,<l>)2-(V$)2)-K($(r/,x)) 



(220) 



where the potential term is as in the text (i.e. it could include a coupling to the curvature 
scalar). Following the standard procedure for obtaining the canonical momentum n(?7,x) to 
$(?], x), and to get at the Hamiltonian density yields: 



Il{r],x) = a^{r])^'{r], x) 

n = ^ + 4^(V<f) + «'(^)m- (221) 

Here conformal time derivatives are denoted by a prime. The generator H of displacements 
in conformal time is the spatial integral of Ti above. Thus the conformal time Liouville 
equation reads: 

^^ = IH,P]. (222) 

If we label the field and its conjugate momemtum in the comoving time frame (i.e. that 
of the text) as $, fl respectively, the results of section 2 are: 



n(f , t) 

H 



a^{t)^{x,t) 
' 112 



d X 



2^3 + i(V'^)'+«'^(*) 



(223) 



The Liouville equation in comoving time can be rewritten in conformal time using the 
relation: d/dt = a{ri)~^ d / drj . After doing this we find that eq. (|115|) becomes: 



dr] 



But 



a{r])H 



d^x 



(224) 



(225) 



Comparing this with eq. (|221[) , we see that we can make the identifications: 11 = 11, \1' = \E'. 
Thus not only is the physics equivalent in both coordinate systems (as must have been the 
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case), but the physics in the two coordinate systems related by a trivial canonical transfor- 
mation. 

We can rewrite all of the comoving time results in terms of conformal time. Some of the 
equations take on a much simpler form in conformal time than in comoving time. This will 
be important when numerical issues are tackled, such as an analysis of the back-reaction 
problem. This work is in progress. 



7 Dynamics of Chiral Simmetry Breaking: Disoriented 
Chiral Condensates 

The proposition has been put forth recently that regions of misaligned chiral vacuum, or 
disoriented chiral condensates (DCC) might form in either ultra-high energy or heavy nuclei 
collisions [^]-[ 56 1 . If so, then this would be a striking probe of the QCD phase transition. It 



might also help explain [p^ , ^ the so-called Centauro and anti-Centauro events observed in 
high-energy cosmic ray experiments [p9| . 



How can we tell, from a theoretical standpoint, whether or not we should expect a DCC 
to form? Clearly, investigating QCD directly is out of the question for now; the technology 
required to compute the evolution of the relevant order parameters directly from QCD is 
still lacking, though we can use lattice calculations for hints about some aspects of the QCD 
phase transition. What we need then is a model that encodes the relevant aspects of QCD 
in a faithful manner, yet is easier to calculate with than QCD itself. 

Wilczek and Rajagopal |^| have argued that the 0(4) linear a-model is such a model. 



It lies within the same static universality class as QCD with two massless quarks, which is a 
fair approximation to the world at temperatures and energies below Aqcd- Thus work done 
on the n = 4 Heisenberg ferromagnet can be used to understand various static quantities 
arising at the chiral phase transition. 

One conclusion from reference ||6^ was that, if as the critical temperature for the chiral 



phase transition was approached from above the system remained in thermal equilibrium, 
then it was very unlikely that a large DCC region, with its concomitant biased pion emission, 
would form. 

The point was that the correlation length ^ = did not get large compared to the 
T^^. A more quantitative criterion involving the comparison of the energy in a correlation 
volume just below with the T = pion mass (so as to find the number of pions in a 
correlation volume) supports the conclusion that as long as the system can equilibrate, no 
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large regions of DCC will form. 

The only option left, if we want to form a DCC, is to insure that the system is far out of 
equilibrium. This can be achieved by quenching the system (although Gavin and M tiller [p^ 
claim that annealing might also work). What this means is the following. Start with the 
system in equilibrium at a temperature above T^. Then suddenly drop the temperature to 
zero. If the rate at which the temperature drops is much faster than the rate at which the 
system can adapt to this change, then the state of the system after the quench is such that 
it is still in the thermal state at the initial temperature. However, the dynamics governing 
the evolution of that initial state is now driven by the T = Hamiltonian. The system will 
then have to relax from the initial state, which is not the ground state of the Hamiltonian 
to the zero temperature ground state. During this time, it is expected that regions in which 
the order parameter is correlated will grow. We can then hope that the correlation regions 
will grow to be large enough to contain a large number of pions inside them. 

The possibility that the chiral phase transition might occur following a quench in heavy- 
ion collisions was explored by Wilczek and Rajagopal in reference |l6l| .They argue there that 
long wavelength fluctuations in the pion fields can develop after the quench occurs. Modes 
with wavenumbers k smaller than some critical wavenumber fccrit will be unstable and regions 
in which the pion field is correlated will grow in spatial extent for a period of time. The 
essence of this mechanism is that the pions are the would-be Goldstone particles of sponta- 
neous chiral symmetry breaking. In the absence of quark masses, the pions would be exactly 
massless when the vf and a fields are in their ground state. However, during the quench, 
the a field is displaced from its zero temperature minimum so that the required cancellation 
between the negative bare mass^ term in the Lagrangian and the mass^ induced through the 
pion interactions with the a condensate does not occur. This then allows some of the pion 
momentum modes to propagate as if they had a negative mass^, thus causing exponential 
growth in these modes. 

Some studies were already reported in which the classical evolution equations of the linear 



sigma model were studied analytically and numerically pT], p3 



In this section we use the techniques developed in the previous section to study the non- 
equilibrium dynamics of the phase transition. In particular we will focus on the formation 
and growth of domains of correlated pions to assess the possibility of large domains leading to 
the production of a large number of correlated pions that may be detected during energetic 
P — P OT heavy ion collisions. 
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7.1 The 0(4) a-Model Out of Equilibrium 

Our strategy is as follows. We will use the techniques developed in the previous section based 
on the functional Schroedinger representation, in which the time evolution of the system is 
represented by the time evolution of its density matrix. 

The next step is to evolve the density matrix in time from this initial state via the 
quantum Liouville equation: 

,n^ = [H,p{t)], (226) 

where H is the Hamiltonian of the system after the quench. Using this density matrix, we 
can, at least in principle, evaluate the equal time correlation function for the pion fields, and 
observe its growth with time. 

Let us now implement this procedure. We start with the sigma model Lagrangian density 



C = ^d^$-d^$-V{a,7f) (227) 
V{a,n) = ^m\t)$ ■ $ + \{$ ■ $y - ha (228) 

where $ is an 0{N + 1) vector, $ = (cr, vf) and tt represents the pions. 

The linear sigma model is a low energy effective theory for an SUi^{2) x S'f/R,(2) (up 
and down quarks) strongly interacting theory. It may be obtained as a Landau- Ginzburg 
effective theory from a Nambu-Jona-Lasinio modelQ]. In fact, Bedaque and Das[^ have 
studied a quench starting from an SUi^{2) x 5'f/R(2) Nambu-Jona-Lasinio model. 

We have parametrized the dynamics of the cooling down process in terms of a time 
dependent mass term. We can use this to describe the phenomenology of either a sudden 
quench where the mass^ changes sign instanteneously or that of a relaxational process in 
which the mass^ changes sign on a time scale determined by the dynamics. In a heavy ion 
collision, we expect this relaxation time scale to be of the order of r ~ 0.5 — Ifm/c. The term 
ha accounts for the explicit breaking of chiral symmetry due to the (small) quark masses. 
We leave arbitrary for now, though at the end we will take = 3. 

Furthermore we will study the dynamics in real Minkowsky time rather than the original 
rapidity variables as first proposed by Bjorken. Recently Cooper et. al.[^ have reported 
on an investigation of the dynamics using the rapidity variables with very similar results. 

Our first order of business is to identify the correct order parameter for the phase tran- 
sition and then to obtain its equation of motion. Let us define the fluctuation field operator 
X(x,t) as 
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a = ^{t) + x{x,t), (229) 
with (j)(t) a c-number field defined by: 



- 1 /,3,Tr(p(t)a(f))_ 

Here Q is the spatial volume we enclose the system in. The fluctuation field xi^y t) is defined 
so that (i) {x{x,t)) = 0, and (ii) x{x,t) = —(t){t). Making use of the Liouville equation for 
the density matrix, we arrive at the following equations: 

m = P(^) = ^l d'x{U^{x)) (231) 

where Ila{x) is the canonical momentum conjugate to cr(x). 

The derivative of the potential in the equation for 7r(t) is to be evaluated at cr = 0(t) + 
x{x,t). These equations can be combined into a single one describing the evolution of the 
order parameter 0(t): 

m + ^ fd'x { ^^S?J!^ ) U=Ht)+x(x,t)= 0- (233) 



J 6a{x) 

To proceed further we have to determine the density matrix. Since the Liouville equation 
is first order in time we need only specify p{t = 0). At this stage we could proceed to a 
perturbative description of the dynamics (in a loop expansion). 

However, as we learned previously in a similar situationfl^, the non-equilibrium 
dynamics of the phase transition cannot be studied within perturbation theory. 

Furthermore, since the quartic coupling of the linear sigma model A must be large (A ~ 
4 — 5 so as to reproduce the value of /^r ~ 95 Mev with a "sigma mass" ~ 600 Mev), the 
linear sigma model is a strongly coupled theory and any type of perturbative expansion will 
clearly be unreliable. Thus, following our previous work[p5|, |39| and the work of Rajagopal 



and Wilczek|^T| and Pisarski ||63|| we invoke a Hartree approximation. 

In the presence of a vacuum expectation value, the Hartree factorization is somewhat 
subtle. We will make a series of assumptions that we feel are quite reasonable but which. 
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of course, may fail to hold under some circumstances and for which we do not have an a 
priori justification. These are the following: i) no cross correlations between the pions and 
the sigma field, and ii) that the two point correlation functions of the pions are diagonal 
in isospin space, where by isospin we now refer to the unbroken 0{N) {N = 3) symmetry 
under which the pions transform as a triplet. These assumptions lead to the following Hartree 
factorization of the non-linear terms in the Hamiltonian: 



6(x^)x^ + constant (234) 

^ 3(x')x (235) 

(vr-vf)^ (2 + ^)(7r2)7f2 + constant (236) 

7?V - ^'{X') + {^')X' (237) 

7?'x ^ (238) 



where by "constant" we mean the operator independent expectation values of the composite 
operators. Although these will be present as operator independent terms in the Hamiltonian, 
they are c-number terms and will not enter in the time evolution of the density matrix. 

It can be checked that when = one obtains the 0{N + 1) invariant Hartree factor- 
ization. 

In this approximation the resulting Hamiltonian is quadratic, with a linear term in x'- 



H^ fr2 



(Vx)^ 



7r 



(239) 



Here H^, H^ are the canonical momenta conjugate to x{.^)i 7f(x) respectively and is 
recognized as the derivative of the Hartree "effective potential" [^, |6^ with respect to (it 



is the derivative of the non-gradient terms of the effective action p5|, |39|). 

In the absence of an explicit symmetry breaking term, the Goldstone theorem requires 
the existence of massless pions, M.t^ = whenever for (p ^ 0. However, this is not the 
case within our approximation scheme as it stands. 

This situation can be easily remedied, however, by noting that the Hartree approximation 
becomes exact in the large N limit. In this hmit, (vf^) ^ 0{N), (x^) ^ 0{1), cjP' ^ 0{N). 
Thus we will approximate further by neglecting the 0{1/N) terms in the formal large N 
limit. This further truncation ensures that the Ward identities are satisfied. We now obtain 



V\t) = <P{t) \m\t) + 4X<f)\t) + 4A(7f2)(t) 



h 



(240) 
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Mlit) = m2(i)+4A02(i)+4A(7f2)(i) 
Mlit) = m\t) + 12X(P\t)+4X{n^){t). 



(241) 

(242) 



The Hamiltonian is now quadratic with time dependent self-consistent masses and Gold- 
stone's Ward identities arc satisfied. 

Since the evolution Hamiltonian is quadratic in this approximation, we propose in the 
Hartree approximation a gaussian density matrix in terms of the Hartree-Fock states. As a 
consequence of our assumption of no cross correlation between x and tt , the density matrix 
factorizes as 

p{t) ^ p^{t) ® p^{t) (243) 

In the Schroedinger representation the density matrix is most easily written down by 
making use of spatial translational invariance to decompose the fluctuation fields x(x, t) and 
TT (x, t) into spatial Fourier modes: 



X{x:t) = -l=^Xkit)^M-^k-x) (244) 
= ;^X^%exp(-ifc-f), (245) 

where we recall that x(x, t) = a{x) — (j){t). We can now use these Fourier modes as the basis 
in which to write the density matrices for the sigma and the pions. We will use the following 
Gaussian ansatze: 



2% ' 2h 



h 



Xk{t)x-k{t) 



h 



Px,,{t) {x-k{t)-X-k{t))\, 



(246) 



-TTfc • 7r_fc H — TTk ■ 7^-k + 



2h 



2h 



h 



(247) 



The assumption of isospin invariance implies that the kernels Aj^^k, B^^^k transform as isospin 
singlets, since these kernels give the two point correlation functions. Furthermore, hermiticity 
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of the density matrix requires that the mixing kernel B be real. The lack of a linear term in 
the pion density matrix will become clear below. 

The Liouville equation is most conveniently solved in the Schroedinger representation, in 
which 

S S 
iiy(x) = —ih— : ni(x) = —ih- — , 

^^^= (^H[U^,li^,x,n-,t] - H[Il^,%,x,^-,t]^ p{t). (248) 

Comparing the terms quadratic, linear and independent of the fields {x ] tt), we obtain 
the following set of differential equations for the coefficients and the expectation value: 



iB 



^X,k 

-Bv k ( 



A 



K,k 



<kit) 



X,k \^X,k 



A 



X,k 



4) + m^it)(f) + 4A0^ + 4A0(7f2(f , t)) - h = 0. 



(249) 

(250) 
(251) 
(252) 

(253) 



^J^T7,k 

AC, A; 



^iA,k 



BjT k ( Aj^ h — A. 



■iT,k 



(254) 

(255) 
(256) 
(257) 



The lack of a linear term in ( 247 ) is a consequence of a lack of a linear term in vr in 
the Hartree Hamiltonian, as the symmetry has been specified to be broken along the sigma 
direction. 

To completely solve for the time evolution, we must specify the initial conditions. We 
will assume that at an initial time {t = 0) the system is in local thermodynamic equilibrium 
at an initial temperature T, which we take to be higher than the critical temperature,Tc ~ 
200 MeV , where we use the phenomenological couplings and masses to obtain T^. 

This assumption thus describes the situation in a high energy collision in which the 
central rapidity region is at a temperature larger than critical, and thus in the symmetric 
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phase, and such that the phase transition occurs via the rapid coohng that occurs when the 
region in the high temperature phase expands along the beam axis. 

The assumption of local thermodynamic equilibrium for the Hartree-Fock states deter- 
mines the initial values of the kernels and the expectation value of the sigma field and its 
canonical momentum: 

A^^k{t = 0) = uj^^kiO) coth[(3huj^^k{0)] (258) 
R ,(t = 0) = -— r^fe^^^^^ (259) 
A^,k{t = 0) = a;^,fe(0) coth['/5acu,,fe(0)] (260) 

B^A' = ^) = --rir^ (261) 
smh[phUn,k[0)\ 

0(t = 0) = 00 ; <P{t = 0) = 0, (262) 

with f3 = l/ksT. We have (arbitrarily) assumed that the expectation value of the canonical 
momentum conjugate to the sigma field is zero in the initial equilibrium ensemble. These 
initial conditions dictate the following ansatze for the real and imaginary parts of the kernels 
A.„^k(t), A^^kit) in terms of complex functions An,k(t) = AR-n,k(t) + iAi-n,k(t) and A^^k(t) = 
AR.^^t) + tAi.^^t)^: 



AR;n,k(t) = AR-n,k(t) coth[j3huJn,k{^)] 

R m = -Ar-ttA^) 
""^''^'^ sinh[/?fic..,,(0)] 

^/;7r,fe(^) = Al-T,^k{t) 



(263) 
(264) 
(265) 



The differential equation for the complex function A can be cast in a more familiar form by 
a change of variables 



A-K^k (^) 



.^^,fc(t) 



(266) 



with obeying the following Schrodinger-like differential equation, and boundary condi- 
tions 



^n,k{t) 
^.,k{t = 0) 







t^7r,fc(0) 



*^,fc(t = 0) = zVcj^,fc(0). 



(267) 
(268) 
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Since in this approximation the dynamics for the pions and sigma fields decouple, we will 
only concentrate on the solution for the pion fields; the effective time dependent frequen- 
cies for the sigma fields are completely determined by the evolution of the pion correlation 
functions. In terms of these functions we finally find 

(7ffc(t) ■ 7r_fc(t)) = ^ I 1 2 coth 

In terms of this two-point correlation function and recognizing that the \I'7r,fc(t) only depends 
on k"^, we obtain the following important correlations : 

(7f2(f,t)) 
(7r(f,t)-7f(0,t)) 

The presence of the temperature dependent function in the above expressions encodes the 
finite temperature correlations of the initial state. The set of equations ( p53|J267| ) with 
the above boundary conditions, completely determine the non-equilibrium dynamics in the 
Hartree-Fock approximation. We will provide a numerical analysis of these equations in the 
next section. 

7.2 Pion production 

In the Hartree approximation, the Hamiltonian is quadratic, and the fields can be expanded 
in terms of creation and annihilation of Hartree-Fock states 

Tfkit) = ^ {Sk^lkit) + Sl.^nAt)) ■ (273) 

The creation aj, and annihilation operators are independent of time in the Heisenberg pic- 
ture and the mode functions "^j^^t) are the solutions to the Hartree equations ( |267| ) which 
are the Heisenberg equations of motion in this approximation. The boundary conditions 
( p68|) correspond to positive frequency particles for t < 0. The creation and annihilation op- 
erators may be written in terms of the Heisenberg fields (|273| ) and their canonical momenta. 
Passing on to the Schroedinger picture at time t = 0, we can relate the Schroedinger picture 
operators at time t to those at time t = via a Bogoliubov transformation: 

Skit) = n,fc(t)4(0) + J-_,fc(t)alfe(0) (274) 



^C^7r,fc(0) 

2kBT 



(269) 



Nh 

Nh 
4^ 



dkk^ I "^nAt) r coth 
sin(A;x 



2kBT 



dkk- 



X 



^^,fc(t) r coth 



^^7r,fc(0) 

2kBT 
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(271) 
(272) 
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with 



1 *7r,fc(t) 


2 " 




2 



1 + 









2 









+ 
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At any time t the expectation value of the number operator for pions (in each fc-mode) 



IS 



Tr 



Tr 



<,(0)-a^,,(0)p(t) 



Trp(O) 

After some straightforward algebra we find 



Trp{t) 



(N^Ait) = (2 I J'+At^Q r -l)(Ar,,fc)(0) + (l T+At,to) 



1 



(276) 



(277) 



The first term represents the "induced" and the second term the "spontaneous" particle 
production. In this approximation, particle production is a consequence of parametric am- 
plification. The Hartree-Fock states are examples of squeezed states, and the density matrix 
is a "squeezed" density matrix. The squeeze parameter (the ratio of the kernels at a time t 
to those at time t = 0) is time dependent and determines the time evolution of the states 
and density matrix. The relation between squeezed states and pion production has been 
advanced by Kogan[70| although not in the context of an initial thermal density matrix. 

Thus far we have established the formalism to study the non-equilibrium evolution during 
the phase transition. A question of interpretation must be clarified before proceeding further. 
Our description, in terms of a statistical density matrix, describes an isospin invariant mixed 
state, and thus does not prefer one isospin direction over another. A real experiment will 
furnish one realization of all the available states mixed in the density matrix. However, if 
the pion correlation functions become long ranged (as a statistical average) it is clear that 
in a particular realization, at least one isospin component is becoming correlated over large 
distances, thus it is in this statistical sense that our results should be understood. 

This concludes our discussion of the formalism we will use to study the non-equilibrium 
evolution of the pion system. We now turn to a numerical analysis of the problem. 



7.3 Numerical Analysis 

The phenomenological set of parameters that define the linear sigma model as an effective 
low energy theory are as follows (we will be somewhat cavalier about the precise value of 
these parameters as we are interested in the more robust features of the pion correlations) 
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M„ ^ 600 MeV ; ^ 95 MeV ; A ^ 4.5 
h ^ (120 MeV f ; 200 MeV . 



(278) 



The above value of the critical temperature differs somewhat from the lattice estimates 
(Tc ^ 150 MeV ), and our definition of A differs by a factor four from that given elsewhere [l6l , 

The first thing to notice is that this is a strongly coupled theory, and unlike our previous 
studies of the dynamics of phase transitions ||16|, ^ we expect the relevant time scales to be 
much shorter than in weakly coupled theories. 

We must also notice that the linear sigma model is an effective low-energy cutoff theory. 
There are two physically important factors that limit the value of the cutoff: i) this effective 
theory neglects the influence of the nucleons, and in the Hartree approximation, the vector 
resonances are missed. These two features imply a cut-off of the order of about 2 GeV; ii) 
the second issue is that of the triviality bound. Assuming that the value of the coupling is 
determined at energies of the order of its very large value implies that the cutoff should 
not be much larger than about 4-5 GeV, since otherwise the theory will be dangerously 
close to the Landau pole. From the technical standpoint this is a more important issue 
since in order to write the renormalized equations of motion we need the ratio between 
bare and renormalized couplings. In the Hartree approximation this is the "wave function 
renormalization constant" for the composite operator vf^. 

Thus we use a cutoff A = 2 GeV . The issue of the cutoff is an important one since 
(vf^) requires renormalization, and in principle we should write down renormalized equations 
of motion. In the limit when the cutoff is taken to infinity the resulting evolution should 
be insensitive to the cutoff. However the chosen cutoff is not very much larger than other 
scales in the problem and the "renormalized" equations will yield solutions that are cutoff 
sensitive. However, this sensitivity will manifest itself on distance scales of the order of 0.1 
fm or smaller, and we are interested in detecting correlations over many fermis. The size 
of the correlated regions and the time scales for their growth will be determined by the 
long wavelength unstable modes|3^ (see below), and thus should be fairly insensitive to the 
momentum scales near the cutoff. The short distance features of the correlation functions, 
such as e.g. the amplitude of the fluctuations will, however, be rather sensitive to the cutoff. 

The most severe ultraviolet divergence in the composite operator vf^ is proportional to A^. 
This divergence is usually handled by a subtraction. We will subtract this term (including 
the temperature factors) in a renormalization of the mass at t = 0. Thus 
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ml{t = 0) + 4A(7r^)(i = 0; T) = m|(i = 0; T) (279) 

where we made explicit the temperature dependence of tt^ and m|j(i = 0;T). Furthermore 
we will parametrize the time dependent mass term as 



mlit) 



Ml 



11 

. c 



exp 



-2- 



- 1 



m + 



Ml 



c 



- 1 



e(-t). 



(280) 



This parametrization incorporates the dynamics of the expansion and cooling processes 
in the plasma in a phenomcnological way. It allows for the system to cool down with an 
effective temperature given by: 



Teff{t) = Texp 



t 



(281) 



where T is the initial value of the temperature in the central rapidity region, and tr is a 
relaxation time. This parametrization also allows us to study a "quench" corresponding to 
the limiting case tj. = 0. 

It is convenient to introduce the natural scale 1 fm ~^ « 200 MeV = Mp and define 
the following dimensionless variables 



0(t) 
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Mpq ; t 
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In these units 
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Thus we have to solve simultaneously the Hartree set of equations: 

+ 4A/^ + 4/AE(0,r) -// 
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(284) 
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E(^,t) = (7r(f,t) ■7f(0,t))/M| 

^^(^) 12-1^^,(0) |2)coth 



47r2 



z 



10T( GeV ) 



(287) 
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with the boundary conditions 

/(O) = /o ; ^ = (288) 
V^,(0) = ^ ; ^(0) = z^W^ (289) 

Finally, once we find the Hartree mode functions, we can compute the total pion density as 
a function of time: 



with (A^7r,g(T)) given by ( |277| ) in terms of the dimensionless variables. 

The mechanism of domain formation and growth is the fast time evolution of the unstable 
modes ||35[|. 

Let us consider first the case of if = 0. Then / = is a fixed point of the evolution 
equation for / and corresponds to cooling down from the symmetric (disoriented) phase in the 
absence of explicit symmetry breaking perturbations. Let us consider the simpler situation 
of a quench (r^, = 0). The equation for the Hartree mode functions (|286|) shows that for 

< 9/2 the corresponding modes are unstable at early times and grow exponentially. 

This growth feeds back on S(t) which begins to grow and tends to overcome the insta- 
bility. As the unstable fluctuations grow, only longer wavelengths remain unstable, until the 
time when 4AS(r) ^ 9/2, at which point no wavelength is unstable. The modes will continue 
to grow however, because the derivatives will be fairly large, but since the instabilities will 
be overcome beyond this time the modes will have an oscillatory behavior. 

We expect then that the fluctuations will grow during the time for which there are 
instabilities. This time scale depends on the value of the coupling; for very small coupling, 
E(r) will have to grow for a long time before 4AS(r) ^9/2 and the instabilities are shut-off. 
On the other hand, for strong coupling this time scale will be rather small, and domains will 
not have much time to grow. 

It is clear that allowing for a non-zero magnetic fleld or /(O) 7^ will help to shut off the 
instabilities at earlier times, thus making the domains even smaller. 

For typical relaxation times (~ 1 — 2 fm/c ), domains will not grow too large either 
because the fast growth of the fluctuations will catch up with the relaxing modes and shut- 
off the instabilities (when Te//(r) < Tc) on short time scales. Thus we expect that for a 
non-zero relaxation time Tr ~ fm/c , domains will not grow too large either because the 
fluctuations will shut-off the instabilities (when Tejj(r) < T^) on shorter time scales (this 
argument will be conflrmed numerically shortly). 
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We conclude from this analysis that the optimal situation for which large DCC regions 
can grow corresponds to a quench from the symmetric phase in the absence of a magnetic 
field. 

Figure (9. a) shows 4AE(r) vs r, figure (9.b) shows T,{z,t — 1;3;5) vs z and and figure 
(9.c) shows Nt^{t) vs r for the following values of the parameters: Tc — 200 MeV ; T — 
250 MeV ; Tr — ; /(O) = ; H — ; A = 4.5 corresponding to a quench from the 
symmetric phase at a temperature slightly above the critical temperature and no magnetic 
field. 

We clearly see that the fluctuations grow to overcome the instability in times ^ 1 fm 
and the domains never get bigger than about ^ 1.5 fm . Figure (l.c) shows that the number 
of pions per cubic fermi is about 0.15 at the initial time (equilibrium value) and grows to 
about 0.2 in times about 1-2 fermis after the quench. This pion density is thus consistent 
with having only a few pions in a pion-size correlation volume. 

The reason the fluctuations grow so quickly and thus shut off the growth of the unstable 
modes so quickly is the strongly coupled nature of the theory. 

The possibility of long range correlations exists if the initial state is in equilibrium at 
the critical temperature. In this situation there arc already long range correlations in the 
initial state that will remain for some time as the temperature factors enhance the con- 
tributions for long wavelength modes since the Boltzmann factor ^ 1/k for long wave- 
length fluctuations. Figures (lO.a-c) show this situation for the values of the parameters 
T ^Tc^ 200 MeV ; = ; /(O) = ; H ^0 ; A = 4.5. Figure (lO.b) shows 
E(2;, r = 1;2) vs z. In this case the number of pions per cubic fermi in the initial state is 
0.12 and reaches a maximum of about 0.17 within times of the order of a fermi/c. The 
pions, however, are correlated over distances of about 4 — 5 fm with a large number of pions 
per correlation volume ~ 50. These large correlation volumes are a consequence of the initial 
long range correlations. This is the situation proposed by Gavin, Gocksch and Pisarski for 
the possibility of formation of large domains, as there is a "massless" particle in the initial 
state. 

We believe that this situation is not very likely as the central rapidity region must remain 
in equilibrium at (or very close) to the critical temperature before the quench occurs. 

To contrast this situation with that of a weakly coupled theory, figures (ll.a-c) show the 
same functions with the following values of the parameters A = 10~^ ; Tc — 200 MeV ; T — 
250 MeV ; Tr — ; /(O) = ; H — 0. Now the fiuctuations are negligible up to 
times of about 4 fm/c, during which time the correlation functions grow (figure 11. b) and 
pions become correlated over distances of the order 4-5 Fermis. As shown in figure (ll.c) 
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the number of pions per cubic fermi becomes enormous, a consequence of a large paramet- 
ric amplification. In this extremely weakly coupled theory, the situation of a quench from 
above the critical temperature to almost zero temperature does produce a large number of 
coherent pions and domains which are much larger than typical pion sizes. This is precisely 
the situation studied previously |^5|] within a different context. 

We have analyzed numerically many different situations in the strongly coupled case 
(A ~ 4 — 5) including the magnetic field and letting the expectation value of the sigma field 
"roll-down" etc, and in all of these cases in which the initial temperature is higher than the 
critical (between 10 — 20% higher) we find the common feature that the time and spatial 
scales of correlations are 1 fm . Thus it seems that within this approach the strongly 
coupled linear sigma model is incapable of generating large domains of correlated pions. 



7.4 Discussions and Conclusions: 

Our study differs in many qualitative and quantitave ways from previous studies. In partic- 
ular we incorporate both quantum and thermal fluctuations and correlations in the initial 
state. In previous studies it was argued that because one is interested in long-wavelength 
fluctuations these may be taken as classical and the classical evolution equations (with cor- 
relations functions replaced by spatial averages) were studied. We think that it is important 
to quantify why and when the long-wavelength fluctuations are classical within the present 
approximation scheme. 

This may be seen from the temperature factors in the Hartree propagators. These are 
typically (incorporating now the appropriate powers of h): 



h coth 



2kBT 



Thus, modes with wavelength k and energies Uk are classical when huk ^ ksT and yield a 
contribution to the propagator 



^coth 



2kBT 



2kBT/uJk 



(notice the cancellation of the h). For long- wavelength components this happens when 

^[rVr,2-i]«i, 



because the "thermal mass" (squared) for the excitations in the heat bath is -f-\T'^ /T'^ — 1]. 
For the phenomenological values of and T^, the "classical" limit is obtained when 
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- 1] « 0.1, (291) 

that is, when the initial state is in equilibrium at a temperature that is extremely close to 
the critical temperature. This is the situation that is shown in figures (3.a-c) where, indeed, 
we obtain very large correlated domains that were already present in the initial state after 
a quench from the critical temperature all the way to zero temperature. 

After an energetic collision it seems rather unlikely that the central region will be so 
close to the critical temperature. If the temperature is higher than critical, in order for the 
system to cool down to the critical temperature (or very near to it) and to remain in local 
thermodynamic equilibrium, very long relaxation times are needed, as the long-wavelength 
modes are typically critically slowed down during a transition. Long relaxation times will 
allow the fiuctuations to shut off the instabilities as they begin to grow and the system will 
lose its long range correlations. 

This was the original argument that discarded an equilibrium situation as a candidate 
for large domains. Furthermore, typical heavy ion collisions or high energy processes will 
not allow long relaxation times (typically of a few fermis/c). Thus we believe that in most 
generic situations, a classical approximation for the long wavelength modes is not reliable in 
the Hartree approximation. 

We should make here a very important point. We are not saying that large coherent 
fiuctuations cannot be treated semiclassically. They can. What we are asserting with the 
above analysis is that within the Hartree approximation, long wavelength excitations cannot 
be treated as classical. The Hartree approximation in the form used by these (and most other) 
authors does not capture correctly the physics of coherent semiclassical non-perturbative 
configurations. 

Thus although the most promising situation, within the model under investigation, is a 
quench from the critical temperature (or very close to it) down to zero temperature, it is our 
impression that this scenario is physically highly unlikely. 

There is another very tantalizing possibility for the formation of large correlated pion 
domains within the linear sigma model and that is via the creation of a critical droplet 
that will complete the phase transition (first order in this case) via the process of thermal 
activation over a "free energy" barrier. The small magnetic field (resulting from the small 
up and down quark masses) introduces a small metastability [ffl|] . The classical equations of 
motion allow for a solution in which the pion field is zero everywhere and a droplet in the 
sigma field (this is the 0(3) symmetric bounce responsible of thermal activation in scalar 



metastable theories in three dimensions O]). Using a spherically symmetric ansatz for a 
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sigma droplet of radius R 

a,iir)^Utanh[M^{r-R)], (292) 

and assuming, for the sake of argument, that the thin-wall approximation is reliable, we 
obtain an approximate form for the energy of the droplet 

E ^ ATcR^f^M^ - —R'hU (293) 

The critical radius is thus (this approximation is clearly reliable only as an order-of-magnitude 
estimate) 

~ 3 — 5 fm . 

Typical free energy barriers are thus of the order of 500-600 MeV. By considering the fluctu- 
ations of the pions around this configuration, it is conceivable (although we cannot provide 
a more convincing argument at this stage) that the unstable mode of the droplet (dilation) 
that makes the droplet grow to complete the phase transition via thermal activation, pro- 
duces a large amount of correlated pions. This scenario, however, requires supercooling (the 
false vacuum to be trapped) which again requires long relaxation times (again unlikely for 
strong coupling). 

As argued above, this possibility cannot be studied via a Hartree approximation which 
only provides a (select) resummation of the perturbative expansion and is probably reliable 
only for short times, before non-perturbative configurations truly develop. 

Thus we conclude that although our analysis provides a negative answer to the question 
of the possibility of large correlated domains near the chiral phase transition, these results 
are valid only within the Hartree approximation of the linear sigma model. There are several 
conceivable possibilities that would have to be studied thoroughly before any conclusions are 
reached: i) perhaps the linear sigma model is not a good candidate for studying the dynamics 
of the chiral phase transition (although it describes the universality class for the static 
properties) ii) there are large coherent field configurations (droplets) that are not captured 
in the Hartree approximation. This possibility is rather likely and is closer to the scenario 
envisaged by Bjorken, Kowalski and Taylor ||57|1 . These semiclassical coherent configurations 
may be responsible for large regions of correlated pions. An important ingredient in this 
latter case must be a deeper understanding of the dynamical (relaxation) time scales, for 
which a deeper understanding of the underlying strongly interacting theory is needed. A 
particularly relevant question is whether such a strongly coupled theory can yield enough 
supercooling so as to produce such a configuration. 
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We believe that these two possibihties must be studied further to give an unequivocal 
answer to the question of large correlated domains. We are currently studying the second 
possibility in more detail. 
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Figure Captions: 

Figure 1: Final contour of evolution, eventually T' ^ oo , T —oo. 

Figure 2(a): i] ys t (notation in the text) for g = 10^^, 7/(0) = 2.3 x 10^^, "^'(0) = 0, 
L = 1. The solid line is the classical evolution, the dashed line is the evolution from the 
one-loop corrected equation (|50|) . 

Figure 2(b): One-loop contribution including the coupling g for the values of the param- 
eters used in Figure 2(a). 

Figure 2(c): Velocity same values and conventions as in figure 2(a). 

Figure 2(d): (gmax(T))^ vs r for the same values as in figure 2(a). 

Figure 3(a): rj vs t with g = 10"'' ?7(0) = 2.27 x 10~^, rj'^O) = 0, L = 1. Same conventions 
for solid and dashed lines as in figure 2(a). 

Figure 3(b): One-loop contribution to the equations of motion including the coupling g 
for the values of the parameters used in Figure 3(a). 

Figure 3(c): Velocity same values and conventions as in figure 3(a). 

Figure 3(d): (gmax(T))^ vs r for the same values as in figure 3(a). 

Figure 4(a): rjvsT with g = 10"^ 7/(0) = 2.258X 10-^ ^^'(O) = 0, L = 1. Same conventions 
for solid and dashed lines as in figure 2(a). 

Figure 4(b): One-loop contribution including the coupling g for the values of the param- 
eters used in Figure 4(a). 

Figure 4(c): Velocity same values and conventions as in figure 4(a). 
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Figure 4(d): {qmax{'T)y vs r for the same values as in figure 4(a). 



Figure 5: Hartree (solid line) and zero order (dashed hne) results for ^^(($^(t)) — 
($2(0))) = V{0, t), for A = 10-12, 1 = 2. 

Figure 6. a: Scaled correlation functions for r = 6, as function of x, V^^^\x,t) (solid 
line), and p(°)(a;,r) (dashed line). A = 10"^^ § = 2- 

Figure 6.b: Scaled correlation functions for r = 8, as function of x, V^^^\x,t) (solid 
line), and V(^\x,t) (dashed line). A = lO'^^, J = 2. 

Figure 6.c: Scaled correlation functions for r = 10, as function of x, T>^^^\x,t) (solid 
line), and V^^\x,r) (dashed line). A = lO'^^^ § = 2. 

Figure 6.d: Scaled correlation functions for r = 12, as function of x, V'^^^\x,t) (solid 
line), and V^^\x,t) (dashed line). A = lO'^^, g = 2. 

Figure 7: Hartree (solid line) and zero order (dashed line) results for ^^((^^(t)) — 
{^^{t))) = 3P(0,r), for A = 0.01, J = 2. 

Figure 8: Scaled correlation functions for r = 4.15, as function of x, T)^^^\x,t) (solid 
line), and V^'^\x,r) (dashed line). A = 0.01, J = 2. 

Figure 9. a: 

4AE(t) vs t. Te = 200 MeV ; T = 250 MeV ; = ; /(O) = ; H ^ ; A = 

4.5 

Figure 9.b: T,{z,t — 1;3;5) vs z for the same values of the parameters as in figure (9. a) 
larger values of time correspond to larger amplitudes at the origin. 

Figure 9.c: N^^^t) vs r for the same values of the parameters as in figure (9. a). A^7r(0) = 
0.15. 
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Figure lO.a: 

4AE(t) vs t. T = Te = 200 MeV ; = ; /(O) = ; H ^ ; A = 4.5 



Figure lO.b: E(2;, r = 1; 2) vs 2; for the same values of the parameters as in figure (10. a) 
larger values of time correspond to larger amplitudes at the origin. 

Figure 10. c: Nt^{t) vs t for the same values of the parameters as in figure (10. a). A^7r(0) = 
0.12. 

Figure 11. a: 

4AE(t) vs t. Te = 200 MeV ; T = 250 MeV ; = ; /(O) = ; H ^ ; A = 
10-6 

Figure ll.b: E(2;, r = 4; 6) vs 2; for the same values of the parameters as in figure (11. a) 
larger values of time correspond to larger amplitudes at the origin. 

Figure ll.c: Nt^{t) vs t for the same values of the parameters as in figure (11. a). A^7r(0) = 
0.15. 
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